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A TECHNIQUE TO DETERMINE THE STRESS 
DISTRIBUTION IN ROTATING ELEMENTS 

ABSTRACT 

A method for measuring In-plane displacement of a rotating 
structure by using two laser speckle photographs Is described. 

From the displacement measurements one can calculate strains and 
stresses due to a centrifugal load. 

This technique Involves making separate speckle photographs 
of a test model. One photograph Is made with the model loaded 
(model is rotating); the second photograph is made with no load on 
the model (model is stationary). A sandwich is constructed from 
the two speckle photographs and data are recovered in a manner 
similar to that used with conventional speckle photography. The 
basic theory, experimental procedures of this method and data 
analysis of a simple rotating specimen are described in the text. 
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I, INTRODUCTION 


Analysis of the stress and vibration characteristics of 
rotating turbine engine components under operating conditions is a 
continuing |.'oblem in the testing and development of today's air- 
plane power plants. J. M. Allen and L. B. Erickson [1] have analyzed 
a free standing gas turbine blade by using the NASTRAN finite 
element code. Correct modeling of the service loads on the structure 
is a major difficulty in such analyses. Other investigators have 
used conventional strain gages to experimentally analyze such 
structures. Strain gage data is usually delivered through slip 
rings, or via radio telemetry, to recording equipment. Problems 
associated with getting these signals to the recorder often plague 
such tests [2]. 

J. P. Sikora and F. T. Menden.v; .1, Jr. observed and analyzed 
vibration-mode patterns of a rotating propeller by means of an 
axially synmetric spinning-hologram technique in 1974 [3]. K. A. 
Stetson also investigated the vibration problem of rotating objects 
by using an image derotator with a hologram interferometry technique 
in 1977 [2], [M- Although holographic interferometry is well 
suited to measurement of nonnal to the surface movement, the 
separation of the in-plane displacements (from the out-of-plane) is 
more complex. Since the meafTement of in-plane displacements is 
usually necessary for strain computation, a simple coherent method. 
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possessing the advantages of holographic interferometry (high 
sensitivity and large field converge of complex shapes) but 
suitable for direct in-plane displacement measurement would be 
desirable. 

Laser speckle photography has emerged in recent years as a 
proficient technique for measuring small in-plane displacements 
[5], [61- This article develops a so called "Sandwich Speckle 
Photography" technique for determining the stress of rotating 
elements. The approach is to take a speckle picture of the 
rotating element with a Q-switched, ruby laser. Next the element 
is stopped in approximately the same position as when the dynamic 
speckle picture was taken and with a new film plate a second 
speckle picture is taken. The photographs are then mechanically 
superimposed, emulsion against emulsion, and carefully aligned to 
produce Young's fringes. This approach supposes that the speckle 
surface will produce the same speckle pattern when Illuminated 
with a pulse laser at different times except for the displaced 
position due to the element rotation. The incentive to develop this 
technique was motivated by the work of Nils Abramson, who developed 
a "Sandwich Holography" technique in 1974 [7]. B. E. Maddux and 
F. 0. Adams developed a Dual -pi ate technique for static speckle 
photography in 1976 [8]. 

In Chapter II, the basic theory of sandwich speckle photography 
is introduced and a mathematical model is developed to explain the 
fringe pattern one observes in the diffraction halo of a sandwich 
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specklegram. The experimental setup and procedures are described 
in Chapter III. Stress analysis of a simple rotating model by 
using this method is presented as an example. 



II. BASIC THEORY AND MATHEMATICAL MODEL 
OF SANDWICH SPECKLEGRAM 

When an ojbect is illuminated with a laser a speckle pattern 
is produced. If the speckle is photographed and then the model 
loaded and the resulting speckle pattern superimposed on the first 
photograph, the double exposed photograph contains infonnation about 
the displacement of the model under load. This information can be 
seen by shining a laser beam through a small area of the negative 
and observing an interference pattern on a screen. The above 
statements are the description of a conventional double exposure 
specklegram [9]. In most of the applications of double exposure 
speckle photography, the specimen will undergo rigid body movement 
in addition to deformation between exposures. In order to cancel 
the effect of rigid body in-plane movement of the specimen, the 
film holder must be translated and/or rotated between exposures by 
a cor»*esponding amount. In most instances the rigid body movement 
is not known in advance of the loading. 

A sandwich specklegram offers the user a method of canceling 
rigid body movement without a prior knowledge of it. Instead of 
employing a double exposure technique, two separate speckle photo- 
graphs are made. The film plates are then sandwiched together and 
data recovery is accomplished as with a conventional specklegram. 
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In-plane rigid body translations or rotations are canceled by 
moving one film plate with respect t-. the other. 

The following procedure is employed to record a sandwich 
specklegram. A glass film plate (typically NAH Agfa 10E75) is 
sandwiched with a clear plate of glass (same size and same thick- 
ness as the film plate) placing the emulsion side towards ihe 
center of the sandwich. A speckle photograph is made of f e un- 
displaced object with the emulsion side of the film plate faced 
away from the surface of the object (the film plate is on the 
lens side of the sandwich). Then, a second speckle photog<^aph is 
made of the displaced object using « new film plate, with an 
identical sandwich procedure except that the sandwich is turned 
around so that the emulsion side of the film plate is towards the 
surface of the object (a clear glass plate is on the lens side o^ 

the sandwich). Both film plates are then chemically processed and 

thoroughly dried. A sandwich specklegram is constructed by placing 
the two speckle photographs emulsion to emulsion such that the 
images are geometrically matched. A description of the mechanical 
positioner to accomplish this task is provided in Chapter III, 

With a conventional double exposure specklegram, in-plane 
displacements result in speckle pairs contained in a plane per- 
pendicular to the optical axis of the camera. This geometry is not 
completely simulated by the sandwi^n specklegram technique. Film 
emulsions are not infinitely thin and glass film plates vary in 
thickness. As a resu't- each speckle of a pair will be recorded 

either forward or behind its' partner. This is true even when the 
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recordings are made such that the sandwichs are put together 
emulsion to emulsion. Therefore, data recovery from a sandwich 
specklegram is complicated by the longitudinal displacement of 
paired speckles due to separate emulsions. A modification of the 
Young's equation derived for ^'onventional double exposure specklegrams 
[9] is required in order to explain c i fringe pattern of a sand- 
wich specklegram and acquire displacement data. 

Consider a point, Pj, contained in an aperture and illuminated 
by light source S as shown in Figure 1. The light wave scattered 
by Pi to a point 0 on a distant screen may be represented by the 
real part of 

El = A Exp [i^il (1) 

where A i ; the amplitude and <j)iisa phase angle. 

l2t a coordinate system XYZ originate at this point. Consider 
a second nearby point P2 located by separation vector B" and of 
coordinates x, y, z. In Figure 1, the optical path 
length of the light f<om 3 to ?i to 0 is denoted by Li and is the 
distance 

Li = 1^1 + 1^1 (2) 

Equation (2) can be written in the following form 

Li = (^ • + (M • (3) 

The optical path length of the light from S to P2 to 0 is denoted 


by L2 and is the distanc 

1-2 “ 1^1 iMl 


( 4 ) 



X 


Screen 
1 / 


?2 (Scattering Point) ^ 



/ 




Figure 1 

Diffraction by Two Points 
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This equation can be written as 

La » (5) 

The change of phase is 

(La “ Lj) (6) 

where x is the wave length of the light source. 

Note that 

^ = (7) 

and 

^ = P^O - B (8) 

Equations (4) through (8) when combined, yield 

^ (P^O • - (PT^.P^S)^ 

- (PlO • P^O)^] (9) 

^ [(pTS - Pis - 2 P^S - B + B - 

* ( PlO • PTO - 2 P^O • B + B • B)*^ 

- (P^.p^S)^- (PlO -PlO)^l (10) 

Because the displacement vector F is very small, relative to PO and 
PS the terms F*B in equation (10) are neglected as higher order 
terms. Then from applying the binomial expansion 

(FI? • p 73 - 2 PI§ •l)*^ « (FTS-FI?)*^- 
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(PiO . PiO - 2 PiO • B)** « (PiO • PiO)** 




( 12 ) 


Therefore equation (10) becomes 



’tis 

(ra-Pi5)‘‘ 


+ 


PiO 


B 


(13) 


The terms in brackets in equation (13) are unit vectors in the 
directions PiS and PiO; therefore, the following change in notation 


is employed 

n. * 




(14) 


PlO 

"o * iPjO • PjO)*® 

So the phase difference in optical path length traveled from S to 
the scattering points and then to 0 is given by 

B- (^3 +Hq) (16) 

Let the direction cosines of n^ be 1, m, n, then 

n^ = -k (17) 

n^ * iT + mj + nF (18) 

The light wave at 0 from the second point, will then be 

E 2 = A Exp [i (4>i + A(|))] (19) 

Consider a small aperature centered in the X - Y olane. For each 
scattering point in the aperature B" = xT + yj light reaching point 
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0 from each point in the aperature will be represented by 

F = X Exp [i(<|»i + \^i)] (20) 

vdiere from equation (16) 

A<fri= ^ (lx + ny) (21) 

The total contribution at 0 from all points in the aperture is 
obtained by integrating over the aperture 

Ui»AExp [i*il // Exp [iA<fril dxdy (22) 

-'-'area 

Consider a second aperture identical to the first but uniformly 
displaced along X, Y and Z by ix, iy and Az respectively, ix, iy 
and Az are assumed sufficiently small so that equation (16) Is valid. 

For each scattering point in the first aperture (at coordinates 
X, y), there corresponds an identical poi nt in the second aoerture '»t 
coordinates x+Ax, y+Ay and Az) emitting a light wave t. ncint 0. 

Each of these waves at 0 may be represented by 

E = A Exp (i(i^i +a+2)1 (23) 


where with equation (16) 


B = (xMx) 1 + (yfAy) j +Azk 


thus 


A <>2 = 


_ -2:r [1(x+Ax) + m(y+A>) + (n-l)Azl 


(24) 

(25) 


The total light distribution at 0 from points in the second aperture 
is 

rr 

(26) 


U 2 *^Exp [i<^ili| Exp [iA* 2 ldxdy 
"^^area 


Comparing Ui with U 2 we find that (Note: Ax, Ay and az are 


constants) 
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U 2 * Ui Exp [-i ^ (lAx+fi6y+(n-l)az)] (27) 

The total wave distribution at 0 from both apertures is 

Uj » Ui + U 2 * F, I l+Exp[-i^lAx+raAy+(n-l)Az)]| (28) 
The intensity observed at point 0 is then 

(29) 

_ ilr 

where Uj is the ccxnplex conjugate of and 
Cl * a proportionality constant. 

If 

= Ui[l+£xp (-18)1 

Then 

» Ui* [1+Exp (19)1 

Fj • Fy = Ui • Fi (l+Exp (-19) + Exp (i9) + 1] 

Ii 

= p- [2+2 cose] 

2 I 1 

= 7T- [1+C0S9] (30) 

where 

Ii = Cl Ui • Ui 

So from equation (28), (29) and (30) we obtain 

I = Ii [1+cos ^ (Ta x+nay+(n-l)Az)] (31) 

Ii is the ■'.itensity due to a single aperture. This intensity is 
noduldted by the bracketed term. Ii typically represents a halo 
imJ the bracketed term represents "Young's fringes" across the halo. 
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Let the X^xiste along the in-plane separation direction of 
the two apertures (then Ay*0). Az represents the "gap" between 
film plates when using the sandwich method. Within the halo the 
observed intensity is Idien 

I oc[l+cos^ (lAx+(n-l)iz)] (32) 

O A 

Fringes occur when 

l&x + (n-l)az= N= + l, 3, 5, — (33) 

Three conditions resulting from equation (33) are considered. 
First, if Az=0 orAx>5^z, then equation (33) reduces to 

l4x = ^ N= + 1, 3, 5. — (34) 

which represents the conventional double exposure specklegrara 

Y 

equation. 1 is usually approximated as ^ where X is the axis 
perpendicular to Young's fringes in the diffraction halo, R is the 
distance from the specklegram to the screen. Second, if Ax=0 then 
equation (33) becomes 

(n-l)AZ=|^ N= + l, 3, 5 — (35) 

Ore can see from equation (35) that the fringe oattern of the 

iecond condition is a set of concentric circles. The spacing 
between these cricles is inversely proportional to az. Third, 
ifAx and AZ are of nearly equal magnitude, the fringe pattern 
for this condition is curved. The curvature of the fringe pattern 
is caused bv the ((n-l)Az] term. 

Data recovery of a sandwich specklegram can be discussed. 
If the curvature of the fringe pattern is very small, then one can 
get the in-plane displacement, ix, from equation (34). If the 
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curvature of the fringe pattern is large. Ax can be evaluated by 
the following procedure. Translate or rotate one speckle photo- 
graph with respect to the other until the fringe pattern becomes 
a set of concentric circles, the longitudinal displacement Az can 
be acquired from equation (35). Then from equation (33) one can 
calculate Ax as 

AX * I - (n-l) Azl N = + l, 3, 5, — (36) 

where 

n=— B 

The direction cosine of the vector from P|to any observation 
point on the screen in the X-direction is 1. The sign of 1 can be 
determined from the diffraction halo finge pattern. If the obser- 
vation points are on the upper half part of the screen, the relative 
values of 1 are positive; if the observation points are on the lower 
half part of the screen, the relative values of 1 are negative. If 
AZ = 0, the fringe pattern is a set of straight lines. From equation 
(34), one knows that for positive values of ax, the fringe numbers of 
the fringes on the upper half of the screen are positive values; and 
the fringe numbers of the fringes on the lower half of the screen 
are negative values. If a x= 0, the fringe pattern is a set of 
concentric circles. From equation (35), one knows that for positive 
value of AZ, the values of the N's are negative; therefore, if ax 
and AZ are positive values, the fringe pattern of the sum of the 
above two fringe patterns is similar to those illustrated in Figure 
2A. Similarly, if Ax is a negative value and a z is a positive value. 
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X 



(A) 


A x>0 


X 



(S) 


A X<0 


Figure 2 

Two Theoretical Fringe Patterns 
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Figure 2B represents the fringe pattern on the screen. If the front 
plate in the sandwich specklegram is a recording of the undisplaced 
object A z will be positive and the relative displacements caused by 
the image recorded on the screen are in a direction pointing toward 
the center of curvature. 



III. EXPERIMENTAL SETUP AND PROCEDURE 


Introduction 

A moving blade was driven by a motor with a rotational speed 
over 3000 rpm. This model was comprised of a thick disk with two 
bars attached on the periphery of the disk. The experiment has 
three main steps. First, a speckle photograph was taken and 
developed of one of the bars of the model while the motor was 
rotating. Let this picture be called the "dynamic picture." 

Second, with the motor stopped, another speckle photograph of the 
same bar was taken, at the same position as that taken in the 
dynamic picture. This picture is called the "static picture." 

Third, these two film plates were sandwiched by a positioner and 
illuminated with a laser; then information pertaining to displace- 
ments and strains in the bar due to acceleration load was obtained. 

Model 

The model, a thick disk with two bars threaded into its 
periphery, was made of steel. The diameter of the disk was 3.50 
inches, and the thickness of the disk was 1 inch. The central 
part of each bar was 3.25 inches long by .50 inches wide by .125 
inches thick. Both ends of each bar were cylindrical in shape and 
threaded with a .50 inch thick diameter and 1.25 inches in length. 
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The end of each bar was threaded into the periphery of the disk 
with a 3.25Hnch by .50 inch plane parallel to the surface of 
the disk. Tfte two bars were 180 degrees apart on the disk. A 
threaded weight (87.52 grains) was attached to the free end of 
each bar to increase the dynamic load and to balance the 
assembly. This model assembly was axially attached to a .75- 
horsepower, 3540 rpm motor. Figures 3A and 38 show the support- 
ing frame, motor and model. The rubber pads under the support- 
ing frame, motor and model were to absorb vibrations. 

Dynamics 

The basic equipment for taking pictures (reference Figure 4) 
consisted of a Q-switched pulsed ruby laser [10], one Spectra- 
Physics, Inc. model 164 Argon Ion Laser with model 265 Exciter, 
two Textronix 549 storage oscilloscopes, one 0-7500 volts DC 
power supply. One pockel cell power supply, one Spectra- Physics 
5 milliwatts Helium-Neon laser, and two photodiode detectors. 

With the motor rotating, one of the bars would intersect 
the argon laser light path, so that oscilloscope I was triggered by 
photodetector I. After oscilloscope I was triggered, t^ msec 
later it triggered oscilloscope II. After oscilloscope II was 
triggered, t^ msec later it triggered the pulse generator; then 
the ruby laser fired. The angle between the position of the bar 
intersecting the argon laser light path and the position of the 
bar illuminated by the ruby laser light was about 180 degrees. 







fsj 

O 


Figure 4 

Schematic of the Equipment Layout 
for Taking Dynamic Picture 
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So if the sum of t, and was very close to the value of the 
half period of the motor, one could get a dynamic speckle photo- 
graph with the bar In the proper position. 

The following were the procedures for setting up the equipment 

1. The ruby laser was aligned [10]. 

2. A black paper target was put about 2 feet In front of 'iie 
ruby lase~, and the laser fired in its random mode [10]; the laser 
left a white spot on the black paper. 

3. A five milliwatts He-Ne laser and two mirrors were placed 
0 .. the top of the ruby laser as shown in Figure 5. The light beam 
which came from the He-Ne laser was parallel to the light beam 
which came from the ruby laser. The two mirrors were adjusted so 
that the He-Ne laser beam hit the white spot on the black paper. 
This aligned the path of the He-Ne laser with the ruby laser. The 
ruby laser is a pulse laser, so it was very difficult to use the 
ruby laser to help one set up the other equipment. But the He-Ne 
laser beam is a continuous beam, therefore one could use the 
He-Ne laser to help set up the other equipment. 

4. The black paper was removed. The motor, model, rubber 
pads, and supporting frames were mounted in front of the ruby 
laser on the floating table. One of the model's bars was 
rotated to the position so that the He-Ne laser beam hit the 
central part of the bar. 



Black Paper 



Figure 5 

He-Ne Laser Setup 
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5. A ground glass was loaded in the film plate holder with 
the ground side of the glass away from the model. A flashlight 
was used to illuminate the model bar for positioning the con> 
vergent lens and film holder. 

6. A divergent lens was positioned so tiiat the He*Ne 
laser light illuminated the whole bar or the part of the bar that 
was to be photographed. 

7. With the He*Ne laser off, the ruby laser was fired. 

Hie divergent light that can« from the ruby laser was centred 
with the divergent light that cmne from the He-Ne laser. When 
the same area of the model bar was illtmiinated, the divergent 
lens was correctly aligned. 

8. The He-Ne laser was turned on again. The image of 
the bar was observed by looking through the back of the ground 
glass toward the convergent lens. A piece of tape was placed 
on each side of the bar image as seen on the ground glass so 
that the image was exactly between these two pieces of tape. 

This helps one decide the delay times of the oscilloscopes 

for taking a dynamic picture. 

9. The mirrors were adjusted so that the argon laser light 
(green continuous light) hit on the central part of the other 
bar of the model. The He-Ne laser was turned off. 

10. Photodetector 1 was in the position where the argon 
laser light could hit directly after the model bar was rotated 


away. 
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The trigger input gate of oscilloscope I was connected to <.he phoio- 
detector I. gate of oscilloscope I was connected to the trigger 
input gate of oscilloscipe II. gate of oscilloscipe II was 
connected to the ruby laser power supply. A* gate of oscilloscope 
II was connected to a pulse generator. Further information regard- 
ing equipment settings is shown in the Appendix A. Photodetector 
I sent a pulse to trigger oscilloscope I, after tj msec delay. 

A^ gate of oscilloscope I triggered oscilloscope ^I; after a t^ 

msec delay, oscilloscipe II sent a pulse frtHn A^ gate to trigger 
the pulse generator and caused the ruby laser to lase. The ruby 
laser included the flash lanq) power supply, pockel cell power 
supply, and photodetect!'** II. The pockel cell was used to inhibit 
lasing until the laser pulse was desired. Photodetector II was 
connected to the input gate of oscilloscope II, so the output of 
the ruby laser could be monitored. 

11. The only remaining problem was how to deteitnine the values 
of time periods tj and t,. The speed of the motor was 3540 rpm. 
so the half period of the motor was 8.47 msec. If one random lases 
the ruby laser, from the output mode shown or* the screen o^ 
oscilloscope II, one knowsthat lasing took place only between .8 
msec and 1.6 msec after the ruby laser's flash lamp initiation. 
Therefore, t^ was set to equal 7.47 msec and t^ was set to equal 
1.0 msec as a first attempt. The steps used in taking 
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a dynamic picture (the steps are discussed later) were then follow- 
ed using a ground glass instead of the regular film plate. When 
lasing took place, looking through the back of the ground glass 
toward the convergent lens, if the bar image was observed between the 
two pieces of tape on the ground glass, the oscilloscopes* time 
delays were correct. If the image didn't come exactly between the 
tspes, the time delay of oscilloscope I was adjust:-! and the 
above procedures were repeated again, until the image was exactly 
between the tapes. 

12. The ground glass was removed from the film plate holder. 

The following were the procedures for taking the dynamic picture: 

1. All equipment was checked. 

2. The motor was turned on and 3 minutes allowed for the 
motor to reach a steady state. 

3. The argon laser was turned on; a black target card was 
put in front of the argon laser so that the black card covered the 
argon laser light. 

4. The oscilloscoce II was reset. 

5. The light in the room was turned off. 

6. The film plate was loaded on the film holder with the 
emulsion side of the plate away from the convergent lens. The name 
and number of the film plate used was Agfa lOE 75 NAH, 

7. The pulse generator and pockel cell power supply were 
turned on. 

8. The ruby laser power supply was turned on and the power 
control switch was adjusted so that the output power reached 1.5 kw. 
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9. The black card was removed. 

10. Reset oscilloscope I and the ruby laser should 
lase. 

11. The argon laser beam was covered again. The pockel cell 
power supply, pulse generator, ruby laser power supply and motor 
were turned off. 

12. The film plate was unloaded from the film holder and it 
was protected fr(Hn the light. 

13. The light in the room was turned on. 

14. One went to the dark room and the film plate was developed. 
Static Picture 

The film holder is shown in Figure 6A and 63. There were 
two small pins fixed on the lower edge of the "window" and one 
small pin fixed on the left-hand side of the "window." There were 
three small screws that protruded from the back of the film holder. 
The lengths of these screws were all the same. When the film plate 
was loaded into the film holder, the face of the film plate 
should tightly rest against the three tips of the small screws 
and the two sides of the film plate should tightly rest against 
the three small pins on the edges of the "window." This pro- 
cedure insured that the film could be loaded in the same 
position each time. 





Figure 6B 

Film Holder (Back Vie' 



29 


After the negative of the dynamic picture was developed, the 
negative was loaded back into the film holder as noted above. 

A light (flashlight) was used to illuminate the model bar and 
its image was viewed through the negative. The bar was slightly 
rotated until the projected image and the negative image coincided. 
The negative was removed, and a clean glass plate with the same 
thickness as that of the negative was loaded into the film holder. 
With the room light off, a new film plate was loaded into the film 
holder with the emulsion side resting against the clean glass plate, 
(the clean glass plate was between the convergent lens and the film 
plate). The ruby laser was fired. The film plate was removed and 
developed. The two speckle photographs (dynamic and static), 
were placed emulsion to emulsion in the positioner (the positioner 
will be discussed in the next section) and illuminated by the 
argon laser light. If a set of fringes was readily observed on 
the screen, then the data collection was successful. If the data 
was not observed, the model bar was not at the correct position. 

By trial and error, one rotated the model bar to the right or to 
the left .004 or .005 inch, and took another static speckle 
photograph and again observed the resulting data, until a set of 
fringes was observed. 
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Positioner 

The positioner was a device for sandwiching the dynamic and 
static speckle photographs. A photograph of the positioner is 

shown in Figure 7. This device had two main parts c a main frame 
and a translation stage. Figure 8 is a photograph of the main 
frame, and Figures 9A and 9B are the photographs of the trans- 
lation stage. The three screws on the main frame could fix the 

position of the translation stage so that the translation stage 
rotates without translation along the screws. Two small 
spring clips control the vertical and horizontal movements of 
the four by five-inch box inside the translation stage. The 
static and the dynamic speckle photographs were sandwiched together, 
emulsion to emulsion, such that the images were almost geomet- 
rically matched. This sandwich of negatives was then placed in 
the positioner such that one negative was on the window of the 
main frame and the other was below the box window of the trans- 
lation stage. The positioner was made of steel. The weight of 

the translation stage made the two negatives press each other 
tightly. The three screws were set at the correct positions so 
that the translation stage was fixed. The translation stage was 
rotated by moving the handle and translated by turning the spring 
clips until the two images were precisely matched. Displacement 
data was observed by passing a narrow laser beam through the 
photographic images. The laser and positioner setup for data 
recovery is shown in Figure 10. 
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Laser and Positioner Setup for Data Recovery 



IV. EXPERIMENTAL RESULTS 

By means of the metJiod discussed in the last chapter, 
displacements can be evaluated and strain information calculated 
for the moving bars. Several of the fringe patterns from the 
collected data are shown in Figure 11. From the pictures of 
Figure 11, the fringes are obviously a set of straight lines. Tnis 
suggest that the weight of the positioner was sufficient to keep 
contact between the two negative emulsions, which is to say that 
Az almost equals zero. 

If a Cartesian coordinate is defined on tne beam (Figure 11) 

the original point is a specific point of the model bar, the X-axis 

is along the center line of tne beam being positive away from the 

center of the disk, and the positive Y-axis is tangent to the disk. 

The coordinates (in inches) of three points on the bar used as 

example points are Point 1-(2.125, .0), Point 2-(2.417, .0) and 

Point 3-(2.75, .0). The magnification factor of the negative 

was 1.5, so the coordinates of the three points on the negative are 

Point 1-(3.1875, .0), Point 2-(3.625, .0) and Point 3-(4.l25, .0). 

To calculate strain at a point, say (x, .0), the following 
procedures were executed. First, the fringe patterns of two 
neighborhood points on the negative: Point A (x+.0625. .0) and Point 

B (X-.0625, .0) were observed at each example point 1, 2, and 3. 

Next, Ax for these two neighborhood points was calculated according 
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-8 

to equation (34) in Chapter Two (where X » 5145x10 cm * 2G2S.59 
-8 

xIO inch). Strain at the particular example point was calculated 
by means of the difference between the^x's of Points A and 8 
divided by .125 inch. 

Calculations for the information shown in Table 1 were made 
using this procedure. The calculations of theoretical results are 
shown in Appendix B. 

It should be noted that the Ax's of Point A and Point B in 
Table 1 are not the exact displacements of Point A and Point B but 
are the sum of the exact displacement and the miss-matched moven»nt 
between two negatives. For each set of data, the miss -matched 
movement of point is the same. The experimental displacenents 

get larger w.i -le observation points are closer to the free end 
of the bar. The . *s of set 1 get larger as the observation points 
are moved toward the free end of the bar. For each point. Ax is 
greater than the experimental displacen»nt; thus, the miss-matched 
movement is a positive value. ITieAx's of set 2 and set 3 get 
smaller as the observation points approach the free end of the bar, 
and one knows miss-matched movenent is a negative value. Because of 
this mismatch the experiirental displacements are not the same as 
the theoretical displacements; but the experimental strains which 
are the slopes of the displacement curves should match the theoret- 
ical strains. 

With a relatively small distance R, one can see from Table 1 
that the difference of fringe spacing between any two adjacent points 



Table 1 

The Theoretical and Experimental Displacement and Strain Data 




Fringe 

Fringe 

Ax of 

A X of 

Experimental 

Theoretical 

Theoretical 

Error of 



Spacing 
of Pt.A 
xlO"iin 

Spacing 
of Pt.B 
xl0‘-Mn 

Pt.A, 

xlO'^in 

Pt. B 
xlO’-in 

Strain 

xlO'^'in/in 

Disg, 
xlO ^in 

Strain 

xlO’^in/in 

Strain 

% 







W 

Set 1 Pt. 
R ^ 

1 

2.4609 

2.5098 

1.8314 

1.7958 

2.8504 

0.7273 

2.9511 

vu 

3.41 

22.25- Pt. 
inch 

2 

2.2656 

2.3047 

1.9893 

1.9555 

2.6973 

0.8395 

2.9251 

7.79 

Pt. 

3 

2.0703 

2.1094 

2.1769 

2.1366 

3.2251 

0.9676 

2.8769 

12.10 

Set 2 Pt. 
R = 

22.25- Pt. 

1 

2.1875 

2.1484 

2.0603 

2.0978 

2.9968 

0.7273 

2.9511 

1.55 

2 

2.3437 

2.3047 

1.9230 

1.9555 

2.6074 

0.8395 

2.9251 

10.86 

inch 

Pt. 

3 

2.5781 

2.5293 

1.7481 

1.7819 

2.6998 

0.9676 

2.8769 

6.16 

Set 3 Pt. 
R- 

1 

8.5937 

8.4375 

2.0595 

2.0976 

3.0510 

0.7273 

2.9511 

3.38 

87.375 Pt. 

2 

8.8281 

8.6719 

2.0048 

2.0409 

2.8898 

0.9251 

2.9251 

1.21 

inch 

Pt. 

3 

9.0625 

8.9062 

1.9529 

1.9872 

2.7410 

0.9676 _ 

2.8769 _ 

4.72 
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Is very small. This is because the image on the screen is small so 
the difference of fringe spacing is not easily evaluated. 

The data of set 3, R - 87.375 inches. The image on the screen 
is much larger and the difference in the fringe spacing between any 
two adjacent points is more readily and accurately determined. 

From Table 1> it can be seen that the experimental strains of set 
3 are closer to the theoretical strains. Thus* the distance from 
negative to screen is an important variable in interpreting the 
data. 

Stress information can be calculated from strain data and 
Hooke's Law. 



V. Conclusion 

The theory and technique for the formation of fringes in a 
sandwich specklegram were developed and extended to the determi- 
nation of surface stress distribution in any rotating element 
problem. This technique can cancel rigid body movements of the 
test specimen so that one can get the deformation of the rotating 
element. The strains are calculated by differentiating displace- 
ment data. Experimental results and theoretical results are in 
good agreement. 

The range of displacement measurement for the conventional 
speckle photography method is from 200x10*® inches to 20,000 x 

C 

10* inches. When the displacements of the specimen are less than 
200x10*® inches, the conventional speckle photography method is 
useless [9]. The disolacement one gets from the sandwich speckle 
photography method are the simi of the actual displacement of the 
specimen and the miss-matched movement of the negatives. Therefore 
if one uses the sandwich speckle photography method, even though 
the displacements of the specimen are less than 200x10*® inches, 
one still can get displacement data and strain can then be calculated. 

When one uses the sandwich specklegram method to analyze the 
stress distribution in any rotating element problem, two things 
should be noted. First, the time delays of the oscilloscopes must 
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be very accurate to obtain a dynamic speckle photograph. Second, 
the distance from the positioner to the screen should be examined 
and optimized for the most accurate analysis of data. 
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APPENDIX A 


EQUIPMENT SETTINGS 


Ruby Laser Bench 

aperture: .109" diameter 

Ruby Lamp Power Supply 

Capacitor voltage: 1.6 kilovolts 

Pulse Generator 

pul se wi dth : 1 . 5 mi croseconds 
mode: pulse delay 

pulse position: 15 microseconds 
output pulse amplitude: maximum 

Pockel Cell Power Supply 

meter reading: 9.5 microamperes 

Oscilloscope I 

mode: B intensified by A 
delay for A: 7.52 milliseconds 
time base A: 5 m-*croseconds/cm 
time base B: 1.0 milliseconds/cm 
sweep setting: single sweep 
coupling B: AC 
slope: positive 
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Oscilloscope II 
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mode: B intensified by A 
delay for A: 1 .0 milliseconds 
time base A: 5 microseconds/cm 
time base B: .2 mi Hi seconds/cm 
sweep setting: single sweep 
coupling B : AC 
slope: positive 
input: 2 volts/cm 
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APPENDIX B 
THEORETICAL RESULTS 

A moving model is driven by a motor with a rotational speed 
3540 rpm. The model, a thick disk with two bars threaded i.ito 
its periphery, is made of steel. The radius of the disk is r. 

The central part of each bar is L long by b wide by t thick 
(referring to Figure 3A, Figure 3B, and Figure 11). Both ends of 
each bar are cylindrical in shape and threaded with a b-inch 
diameter and Lt in length. One end of each bar is threaded into 
the periphery of the disk. The two bars are 180 degrees apart 
on the disk. A threaded weight M is attached to the free 
end of each bar to increase the dynamic load and to balance 
the assanbly. 

If the model and the threaded weight are assumed tobe homogeneous 
and mass distribution is uniform, then the force acting on the 

bar is 

Fx =M( R + L + Ji L^) W2 + p . L. (lil ) (R + L + L<-)W2 

t ^4 

+ b ' t ML - X) ' p ’[ R + X + Sg (L - X) ] W2 (Al) 

Fx = (M + ^ ir p Lt b2) (R + L + Jg Lt) 

4 ^ 
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+ ^ b t (L - X) (R + 4 L + ^ X) W- (A2) 

where 

W is the angular velocity of the rotating model 
p is the density of the model 

. fx_ * 

A bt (A3) 

where 

is the normal stress of the bar 
A is the cross section area of the bar 

According to Hooke's Law one can get the strain and displace- 
ment by 


0 



and 

U * c d X 

\ X 


(AS) 
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ABSTRACT 


Laser speckle interferctnetry Is used as an experimental 
technique to measure in.plane surface displacsnent components of a 
deformed solid. A general theory of laser speckle interferometry 
which includes both the contributions of in.plane and out of plane 
deformations is developed and verified through several experimental 
examples. The basic information consists of a double exposure 
photographic record of a sheared wavefront of the image of the 
surface points of a body in a reference and deformed configurations. 

Data analysis in speckle interferaretry yields a map of the 
surface displacement components along a specified closed contour 
for the region of interest. The experimental displacement values 
are then used as input data to numerically calculate, with a 
high degree of resolution, the state of stress at any desired 

interior point using the boundary- integral equation method. This 
method is based on the use of an integral constraint equation 
which relates boundary displacmnents and corresponding boundary 
tractions in a well -posed boundary value problem. 

Several examples are presented to illustrate the application 
of the coupled laser speckle interferometry and boundary- integral 
solution technique to two-dimensional elasticity problems. This 
coupled technique between laser speckle interferometry and the 
numerical boundary- integral equation method proved to be an 
accurate and efficient procedure for solving significant engineering 
problems. 
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I. INTRODUCTION 

Recent advances in the area of experimental stress analysis 
utilizing coherent optics have been through th*» development of holo- 
graphic and laser speckle interferometr 7 . These opvical techniques 
"have been applied successfully to the measurement of the surface 
displacements of deformable solidsCl,2,3,4,5,6]*. While holographic 
and speckle interferometry measure surface displacements, speckle 
shearing interferometry determines directly the derivatives of 
surface displacemei> s [7,8]. 

Original applications in holography suggested that the modern 
techniques in using the laser vould possibly change the method of 
measurement of strains. However, this has not been the case, 
although many important engineering examples have been demonstrated. 
Difficulties exist in holography because of the sensitivity of the 
measurements which require vibration isolation and the ability to 
seoarate displacement components from a single hologram. Intro- 
duction of the laser speckle effect has provided for a technique of 
surface displacement measurements without the need of vibration 
isolation. In addition, this technique is a direct measure of in- 
plane displacement components, thus the separation of components is 
accomplished easily. All coherent optic?"' methods are a measure of 
surface displacement or its derivatives on the boundary of a solid, 


* Numbers in brackets refer to the list of references. 
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thus the need still exists to calculate strain and stress 
components at any location fran the surface displacement data. 

The boundary* integral equation (BIE) nsthod is well 
suited to the analysis of laser speckle surface displacement 
data. This method derives its foundations from the application 
of classical potential theory [9], and was used to provide 
numerical solution for the two dim onal Laplace equation 
by Jaswon [10] and Symm [11]. The numerical BIE method has been 
applied to two-dimensional problems of elasticity by Rizzo [12] and 
subsequently extended to three-dimensional problans by Cruse [13, 
14,15]. In addition, boundary- integral methods have been used 
successfully to solve problems in el astodynami cs [16,17], elastic 
inclusions [18], elastoplastioity [19,20], fracture mechanics 
[21,22], thermoelasticity [23], anisotropic elasticity [24], time- 
dependent inelastic deformation [25], axisynmetric problems [26], 
and plate [27] and shell problems [28]. Also, the possibility of 
coupling the BIE method with the finite element method has been 
snown [29]. 

Basically, the BIE method for elasticity combines a reciprocal 
work theorem and a suitable singular soltuion to the governing 
di^erential equation to construct a boundary formula which is a 
relation between boundary displacements and corresponding boundary 
tractions i'^ a well-posed prcol-* Since either of tnese boundary 
quantities, in principle, deternnes the other, the formula provides 
a constraint between them which generates a itt of simultaneous 
integral equations from suitable boundary data. The unknowns in 
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the equations are boundary tractions or displacements directly. 

An outstanding conceptual advantage of the boundary formula 
approach is that it applies to three-dimensions as well as two 
with similarly defined integral equations arising in each case. 

The key feature of the boundary integral equation is that it 
couples boundary data for the given problem. As a result, the 
problem size is reduced by one dinension (with a corresponding 
decrease in computer core requirements). Furthermore, since the 
nunerical discretization is performed at the boundary and not 
throughout the interior, a high resolution for the interior solu- 
tion may be obtained at any desired location. This is particularly 
crucial for analysis of problems with significant stress gradients. 

In principle, the optically measured surface displacements 
beccxne input data needed for the BIE to calculate the unKnown 
tractions at specified points on the boundary. In regions of 
high stress g’^adients or where loading is unknown, a closed 
contour (two-dimensional problems) may be defined around the region 
of interest and the displacement comporants are then measured 
experimentally at specified nodal points on that contour. Stresses 
at any interior point in that region may be obtained accurately 
by comparatively simple quadrature once all the contour tractions 
and displacements are determined. 

The coupling of these two solution techniques becomes more 
compatible wher, the displacement data can be stored in the computer 
directly without the need of time consuming data analysis. Auto- 
mated data analysis systems have been developed [5] for use in 
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laser speckle interferometry which will record and store in memory 
in a computer the displacement components around a specified closed 
contour. 

The subject of this dissertation is the extension of the use 
of the numerical boundary integral technique for elasticity problems 
through the coupling of laser speckle interferometry measurements. 

A general theory of laser speckle interferometry which includes 
both the contributions of in.plane and out of plane deformations is 
developed and discussed in Chapter II. The mathenatical formulation 
of the BIE method for elastostatics and the numerical solution 
for plane problems are discussed in Chapter III. A computer aided 
data reduction systan is also discussed, and several example 
problems are presented to illustrate the coupling of the two 
solution techniques. 
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II. GENERAL SHEARING SPECKLE INTERFEROMETRY 
2.1 Introduction 

When a diffuse surface of a structure is illuminated by coherent 
light, a grainy speckle effect is imaged by the eye or film plane of 
a canera due to the randan interference of light reflected from the 
surface. This phencn^non is known as the spedcle effect. Laser 
speckle interferometry is a relatively new experimental technique 
which shows pranise of elimination to many of the vibration prdJlems 
experienced thus far in holographic interferometry. The basic method 
utilizes high resolution photographies of a surface which is illumi- 
nated with coherent light (laser). The nwthod of recording displace- 
ment information is to photograph t jrface in both a reference and 
deformed configurations which results in a permanently stored whole 
field record through interference fringes of a deformed surface. 
Boundary point in-plane displacement measurements are obtained by 
speckle interferometry and derivatives of the boundary displacements 
are obtained through shearing interferrretry. 

Measurements of in plane displacements and the out c .’.ne deriv- 
atives are determined by the data analysis of a double exposed photo- 
graphic record. Values of displacements and derivatives are determined 
through point-wise and full field Fourier filtering of the film. Point 
wise filtering is a method of analysis in which the laser beam is 
passed through a point on the film plane. As the beam passes through 
the film the defomeo ard undeformed speckles produce a complex dif- 
fraction pattern which causes an interference effect on a viewing 



6 


screen. A diffraction halo modulated by light and dark bars of light 
Is produced where the distance between the bars is Inversely prt^or- 
tional to the distance between the undefom«d and deformed speckle 
on the film plane. A normal to the light and dark bar pattern Indi- 
cates the direction of displacement of the speckle. Displacements 
are thus obtained at a point and each point on the film Is analyzed. 
Whole field Fourier filtering Is used to obtain the derivatives. In 
this form of data analysis, an interference fringe Is the locus of 
points of constant values of the derivative. Thus, in general, a 
photographic record of a double exposed laser speckle interferogram 
contains both displacanent and derivative information. Each method of 
data analysis will be discussed and illustrated by an engineering 
example. Also, the results imply the utilization of a computer based 
data analysis system which will also be discussed In Chapter IV. 

2.2 General Theory of Frinoe Formation 

This section presents the general theory of shearing specxle 
interferanetry. A scharatic diagram of the experimental apparatus for 
shearing speckle interferometry is shown in Figure 2-1. The optical 
wedge in the system shifts the focused image in the direction of the 
shear. Thus, a point p is imaged at (x,y) and (x+6x,y+oy) in the film 
plane, and a neighboring point p^ is also imaged at (x,y) in the film 
plane. In order to obtain fringe data, a double exposure technique 
is used, where a photograph is taken of the body in some reference 
configuration. After defomation, another exposure is recorded on the 
same i Im. 

With th-» body in the reference configuration, the light 





Fifiure 2-1. Shearinti Speckle Intcrfernnetrv 
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amplitudes at the film plane for both points p and Pj^ can be described 
in the following mar*ner 


Sp * i^(x^) expCi9p(x,y)] 

* L (x,y) exp[ia. (x.y)] 


(2-1) 


where a (x,y) is the amplitude modulus and e (x,y) is the phase and 
both are, in general, functiOTS of the coordinates of the film plane. 
The total intensity for the first exposure can be written as 
I, * (^ * I ).(£* + I* ) 

where E* and E* denote the complex conjugate of the light amplitudes, 
or 

Ij * ap(x,y) + 3p^(x,y) + 2ap(x,y) ap^(x,y) cos ♦(xor) (2-2) 

where 

•©■(x.y) « e_(x,y) - 9 (x,y) 

H 

If the body is deformed and the points p and Pj^ are displaced to 
points p' and pj^, then the image of point ? on the film plane will be 

shifted to new coordinates x* and y'. The light amplitudes in the de- 
formed configuration can be written in the following form 


Ep. = ip(x',y') exp[i{9p(x' ,y’) + a9p}] 

'5 (x',y’) exp[i{9 (x',y‘) + ae >] 
Pg ?2 ?2 P2 


( 2 - 3 ) 


i9p and i9p account for the phase change in each field due to motion 
of the object relative to the illumination field. The total intensity 
for the second exposure is 

*e;,) 
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or 

Ig.® ap^(x' ,y') + ap^^(x'.y') 

+ 2 ap(x*,y') ap^(x' ,y')cos['j'(x' ,y‘) * (2-4) 

where 

t(x'.y) • 9p(x'.y‘) - 9p^(x.^.) 

and 


A A0_ - A6_ 

P Pg 

So the total intensity expression for both exposures is 

^TOTAL " h* h‘ 2 3^(x.y)ap^(x.y)cos-»(x,y) 

+ ap*(x',y') + ap^(x',y‘) 

+ 2 a (x',y') a (x',y') cos['?(x' ,y' ) + a«-] (2-5) 

K K2 

Let the displacement vector be written in the Following 
manner (see Figure 2-2). 


Dp'sU = ue + ve + we 
'^p “p X p'y Vz 

then the displacement vector in the film plane can be written as 

[p> (x.y)]FILM = ' “p, ®x * ’Pf Vx * ’p^y' 

where M * magnification of camera, typically M<1 and Up = MUp, 
Pf P 


The coordinates at the film plane for the second exposure are x' * 

X + u^ and y' = y + v^. Now the expression for the total intensi ty. 
Equation (2-5), can be written in the following form 
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ItotAL * * 2 ap(x,y)ap^(x,y)cos-e-(x,y) 

♦ 2»p* > *P2 ' **“Pf V ' ’‘"“Pf 

* A (2«6) 

Interference fringes are obtained by taking optically the Fourier 

transform of the amplitude transmission function of the transparency. 
This can be done by several rethods depending on the data analysis. 

The amplitude transmission function g(x»y) is linear for the ranges 
of interest and can be expressed as 

g(x.y) • a b (2-7) 


where a and be are constants. 

Substituting the expression for from Equation (2-6) into Equa- 

tion (2-7) an expression for the anplitude transmission function is 
obtained as 

g(x.y) * a + b {a_^(x,y) + a_^(x^^) + 2 a_(x»y)a_ U.y)cos «.(x,y) 


'1 


'1 


+ 2 a_(x+u_ ^*v )a. (x+u_ ,y*^v^ )cos['?(x+u^ ,y+v^ 

P P2 Pf Pf Pf Pf 

+ i «•]} (2-8) 

The procedure developed by [ 6 ] is used to express the light 
amplitude 6 (uj,o) 2 ) in the Fourier transform plane (Figure 2-3) which 
is proportional to the Fourier transform of g(x,y) times a quadratic 
phase factor, i.e. 
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Parallel 

Light 

Beam 


Figure 2-3. Arrangement For Optical Fourier Transform 
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• expC — 5^-2-] F[g(x,y)] (2-9) 

P Cg(x,y)3 denotes the Fourier transform of the function g(x,y), 
rather 

W 

p Cg(x,y)] « II g(x,y) exp[-i(u),^x + t^j^y)] dxdy 

and d is the distance between the transparency and the screen. In 
the transform plane, the vectorial spatial -frequency variable is de 
finded as 


ui ® u, e + u>« e ; 

i As t. JTj 


kx 

s 

“1 


(2-10) 



The coordinates in the transform plane are denoted as (x^,y^) and the 
••nstant k « where x is the wavelength of the light. 

Substituting g(x,y) from Equation (2-7) into Equation (2-9) 
yields 


id(u;?+u«) 

S(u 2 »u 2 ^ ® exp[ ■] F [a + b 

, 2 

id(uj>t.)«; 

“ exp[ . 1 2 ~i. 1 . - 1 . -] {a6(u)^ ,(i» 2 ) * b 
where o((*)^,u 2 ) is tne delta function defined as 

CO 

5(id^,.i)2) = [[ exp[-i(u^x+u) 2 y)] dxdy 
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The delta function represents the (idealized) point focus of the 
illuminating bean. It contributes to the amplitude G(u)j,ij 2 ) only at 
the center (uj»u 2 » 0 ); in practice it is a small area around this 
point. Outside this small area the intensity lG(u)j,u 2 )r in the 
transform plane is 

2 

I(uj^,u2) * !G(w2,o) 2)| * G(a)j,u2) G*(uj,u2) 

“ ^ I (2-11) 


Using the identity 

cos('?+a«‘) • cos « cos sin 'y sin 

in the expression of Equation (2-6), the Fourier transform of 

ItOTAL written as 

F [ItotaJ ’ Ap(u)*Ap(u) a + exp [i(;^. Up^)]} + Ap^(u)*Ap^(t;) 

A (ui)„A ((j) exp[i(u. U )] + 2A (u))*A ((i))*C(ui) 

+ 2A (u)*A (u)*C(u) exp[i('j. U )]cos 


- 2A_(u)*A (u)*S(u) exp[i(u. U_ )]sin a «• (2-12) 

P P2 * 

F[ ] indicates the two-dimensional Fourier transform with re- 
spect to the variable u; the capital letters A., A_ and A_ are the 


P Pj K2 

Fourier transforms of the functions a_, a and a respectively; 


P 

denotes the convolution integral ^such as) 


P2 


A_(u)*A (u)*C(w) » F [a.(x,y)a„ (x.y)cos f(x,y)] 


'1 


P-. 


and 


C(ti) = F[cos '?(x,y)] 
S{Z) » F[sin '?(x,y)] 
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In writing Equation (2-12), shi^'t theorem has been used (such as) 

F [a-^Cx+u ,y+v )] * exp[i(u.U )] F[a ^(x,y)] 

P Pf Pf Pf p 

and the convolution theorem to describe the transfom of the various 
products in Equation (2-6). Also Up has been assumed to be a slow 
varying function of x and y. 

Several experimental examples will be discussed to illustrate 
the validity and the data analysis of Equation (2-12). 

2.3 Experimental Verification of the General Shearing Speckle 
Interferometry Equations 

The experimental problem chosen to verify the results of the 
general theory is the flexure of a rectangular beam by a terminal 
end load. Let the beam have the end x = 1 fixed and the end x = 0 
have a distribution of forces statically equivalent to a single force 
Wy directed along the negative y-axis as shown in Figure 2-4, Take 
the X-axis along the center line of the beam, and y-axis as orthogonal 



Figure 2-4. Rectangular Plexiglas Beam 
with a Transverse End Load 
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intersecting at the centroid of the free end. Material fcr tne beam 

5 

was Plexiglas with E ® 4.56x10 psi and » 0.37. The transverse end 
load was governed by the displacement of the free end 6 because of the 
constraints on the pointwise filtering of the experimental data. For 
completeness of the discussion the theoretical displacement components 
from the stress function solution for this problem will be listed here. 


2 

u = C^x^y - ^v+2)y^ + {c^(l+v) - Y 
v,^C.|xy2. Ix^.ll^x- ll3] 


(2-13) 

(2-14) 


Example 1 - In-plane Displacemer. . 

The geometrical configuration for in-plane ir.e^s_; . ents is shown 

in Figure 2-5. A point P on the surface of the b-. . . is illuminattu- 

by a laser located at position S. The surface is imaged by a camera 

located at position 0 with film plane coordinates parallel tc “*'e 

x,y plane. For this example, the optical wedge is removed, thus, the 

lateral shear is zero, which yields, A = A = 0 in Equation (2-12). 

Pi ?2 

Furthermore, since only a small area is illuminated in the film plane 
(see Fgirue 2-6), the displacement in the neighborehood is assured 
constant, the total intensity. Equation (2-11), reduces to the 
following form 

I (u 2 ^,oj 2 ) = 2b^jAp (d)) ^ Ap (u) 1 ^[l+cos(u . Up (2-15) 


Equation (2-15) agrees with the simple case considered in Reference 
[ 53 . 
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Figure 2-5. Geometrical Configuration For The 
First Exoerimental Example 
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Only terms that relate to tJw interference effects in Equation 
(2-15) will be considered in the data analysis. Interference fringes 
will be defined when I * 0 or [l+cos tl.Up ] * 0. Fringes 

are determined by the technique of pointwise filtering as illustrated 
in Figure 2-6. For the small region of illumination of the trans- 
parency in the neighborhood of P, the displacenent is assumed constant 
and fringes are formed when 


u>.U„ = (2n-l)ir 


n • 1,2,3, 


For most problens of interest only the first fringe is observed, 
therefore , 



Displacement components are obtained from the following expression 

2ir / ^S ^S \ » 

T % * T • 'V* * ’"A’ ’ ’ 

Usually 


* ®x* ®y^ “ ®y* 


therefore 


X* Un y* V 


\d 


P, " ^ % T 


(2-16) 


The schematic for determining the components u and v from 

Pf 

Equation (2-16) is illustrated in Figure 2-7 by measuring the 
distances x^ and y^ between the center and points of intersection of 
the first fringe with the x and y axes respectively. 




Figure 2 - 6 . Data Analysis for Point-Wise Filtering 
in ln-p»ane Measurements 


CO 
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Central 

Bright 

Spot 


Figure 2-7. Schematic for Determining the 

In-plane Displacements Components 


The displacement ccmponents for the film plane are then given by 



and 



Figure 2-8 is a graph of the results of in-plane measurements for 
this experimental example. As can be observed, the experimental 
results agree very well with the theory. In addition, this method 
of data analysis is amenable to autanated data analysis as wi**! be 
described later in Chapter IV. 


Example 2 - Shearing Speckle .Measurements of 
Displacement Derivatives 


The experimental arrangement for the example problem used to 
illustrate the shearing speckle measurements is shown in Figure 1-9. 
An optical wedge is placed in front of the camera at 0. ;iso for 



Dimensionless Distance (~p) 

Figure 2-8. Comparison of Theoretical and Experimental Values 
of Oeani Deflection Alonn the x-axis 
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Light 

Sourc’ 


Figure 2-9. Geometrical Configuration For The 
Second Experi frontal Exairole 
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this example u„ ® = 0 while 0 and » A„ = A„ = A since 

P P P P Pj^ P2 

shearing spcdcle interfercjmetry is a phase phencmenon. Equation 
(2-12) reduces to 

F ^TOTAl] * ^ A(iu)^A(<jj) + 2A(u»)*A{u)^^C(<u) (1 + COS 

- 2 HZ)J^C)^S{Z) sin (2-17) 

Again Equation (2-17) agrees with Reference [5]. 

Derivative information is related to the minimum of I((u^, <*> 2 ), 
Equation (2-11), which will be minimum - upon substitution of 
Equation (2-17)-when (I + cos a-») and sin a-» are both zero. There- 
fore, A'»‘ » (2n-l)ir will define a fringe. The phase change A-a- can 
e written in the following form [l] 


~ C^s + ^ 0 ] . (P2?2‘ - PP') 


(2-18) 


where ps 




m. 


and po = 


are unit vectors 




representing the directions from point P 0 ^ the object to the light 
source and camera respectively. Further, as shown in Reference [5], 
Equation (2-18) can be written in the following form 


* ("s * "'o’ 'I 7 * It 

Mn, * Ho) (|2 ox * |2 oy . II ,Z)] (2-19) 


Fringe order can be determined from the following expression 
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♦ ("s * "o> <17 “ * I 7 ll ^»i] <2-201 

Fringe order data for the shearing speckle Is obtained by 
taking optically the Fourier transform of the double exposed trans- 
parency as shown in Figure 2-3. Data for this example is the 
following 

^0 • “"o ’ "0 ’ ^ 

and 

Ig * -sins , a 0 » f'j * coss where s*19.59“ 

For a shear of value Ax * 0.0899 in. in the negative x direction 
only and assuming that the only non-zero displacement canponent is 
Wp, then the fringe order n is expressed as 

n a P sins I7 + (1 + cosS) |fj x ^^-21) 

where \ is the wavelength of the laser light used * 6328 A“. 
Comparison of the theoretical and experimental results .for the case 
of free end deflection of 0.005 inches are shown in Figure 2-10. 

As in the in-plane example the results are in good agreement. 

Both experimental p^amples illustrate the combined effects of 
displacement measuremer.wS in laser speckle interferometry. Also the 
measurements are not constrained to vibration isolation environments 
as in holography, therefore, measurements on prototype systems are 
possible. Possibly one of the most important advantages of the 
soeckle technique is the ability to automate the data analysis by a 
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computer based system, in particular, the in-plane measurements. 
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III. THE BOUNDARY- INTEGRAL EQUATION 
METHOD FOR ELASTICITY PROBLB<1S 

3.1 Brief Introduction 

The boundary-integral equation (BIE) metiod, with its origin in 
classical elasticity, has only in recent years been developed and 
applied to solve significant problems in solid mechanics. The 
initial applications of the BIE metb*'d to solid mechanics problems 
have taken place within the past twelve years [12]. Since then the 
potential of this powerful technique as a basic tool for solving 
engineering prcblems has been recognized and the number of problem 
areas in solid mechanics that have been solved using this numerical 
technique is large [30]. 

The BIE technique is based on the use of an integral con- 
stra-nt equation which relates boundary displacements and 
corresponding boundary tractions. One advantage of this method is 
that numencal discretization occurs only on the boundary of the 
body being anc.’yzed, which reduces the dimension of the problem by 
on?. No approxi;naticn to the field equations are necessary as all 
approximations are maoe at the boundary, hence a high resolution 
for the inte*”ior solution may be dbtained at any desired location. 
This is particularly important for the analysis of problems with 
significant stress risers. The BIE method is also independe“. of 
boundary shape and connectivity of the body. 

Following this section is an extensive review and a thorough 
discussion of the mathematical formulation of the BIE method as 
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applied to elastostatics. Formulation of the numerical solution 
tedinique for plane problans is discussed in Section 3.3. 

3.2 Mathematical Formulation of the BIE Method 

A comprehensive review of the mathematical basis and formula- 
tion of the BIE technioue for two and three dimensional elasticity 
is essential for conscious and better understanding of this powerful 
technique and will be documented in this thesis for completeness of 
the discussion. 

This discussion is restricted to the analysis of classical 
elastostatic problems for which the material may be taken as 
isotropic and homogenous. The notation used is the usual cartesian 
tensor notation with implied summation on repeated indices and 
partial differentiation denoted by the comma-index. 

Equilibrium equations can be written as, 

Cijj + b. = 0 (ij * 1.2,3) (3-1) 

where b^ is the body force. The stress components and displace- 
ment gradients are related by Hooke's Law 


<3 



»{u 


1 .J 


u . .) 

J.V 


(3-2) 


where u is the shear modulus of the material and v is Poisson's 
ratio. Substitution of Equation (3-2) into Equation (3-1) results 
in the Navier equations of eoui libriam in terms of the displacements 






(3-3) 


The solution of this differential equation must satisfy appropriate 
boundary conditions. The displacement boundary value problem 
assumes knowledge of the displacement on the entire surface S 
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u^(x) • ij^(x) (3-A) 

where x denotes the orthogonal cartesian coordinates Xj^,X 2 »X 2 . The 
are the prescribed values of surface displacements. The trac- 
tion boundary value problem is 

ti(x) * o^j-n^ ■ t^(x) (3-5) 

where the vector function t^(x) is prescribed for x^S. The unit 
vector n. is the outward normal vector for the body R. 

J 

The surface tractions and disp'acements associated with any 
two solutions to the governing equations of elasticity in a region 
R can be related through Betti's theorem of elastic reciprocity [31]. 
Betti's second theorem states that for continuous, finite stresses 
and zero body force 

tju^dS = t^ujdS (3-6) 

V J 

S S 

where S is the bounding surface of R. 

If body forces are included in the stress field Betti's second 
theorem can easily be shown to be of the form 

j tju^dS + x|u^dV = ' -.fu|dS + | X^ujdV • (3-7) 

S R S R 

1 2 

for body force fields Xt, which are finite and piecewise 
continuous in R. 

The development of the boundary integ j. approach proceeds by 

2 2 

setting cne of the stress states u^ and t^ in Equation (3-6) jqual 
to the displacements u.* and tractions t^* corresponding to the 
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fundamental solution to Navier's equations, ujau. and tjst^ are 
displacements and tractions for the other stress state. It has 
been shcjwn in Reference [32] that the fundamental solution to 
Navier's equations corresoonds to Kelvin's problem. By letting the 
distance between the field point q with coordinates and y^ 

and load point p with coordinates Xp and x^ be given by 

r = C(y^. - X.) (y^ - x^.)]^^^ (3-3) 

the well known solution [31] to Kelvin's problem of the point load 
of unit magnitude in each of the x^. directions in an infinite body 
is 


i*''’'"' ' r * '-.I'-dlej (3-9) 

where e.'s are a set of unit vectors in the x. directions. In 

J w 

Equation (3-9) and in what follows all differentiations are with 
respect to the field point q, i.e.. 


3r 1 / 

*^»i ° 3y^ ” r ^-^i ” ^i' ax.. 




I?=l7, "i ’r - Xi)"i 

where the normal is evaluated at o also. The tractions on an 
arbitrary su*'face with outward normal n. around the point p 

J 

corresponding to the fundamental solution are given by 

-n.r . + n.r .1 e. 

J »i »J I J 


1 


(3-10) 

(3-11) 


/ 1 « N 
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From Equations (3-9) and (3-12) second-order tensors can be 


found 


(3-13) 


(3-14) 


(3-15) 


* 

u^ (p.o) » Uj^(p,q) Oj . 
t^ (Ptd) » T^^(p.q) 6j 
wh<*r [13] 

■ TSSrtW r ‘ij ♦ ^1^j] 

and 

- n.r . + n.r .1 

The first index in Uj^.(p,q) and T^^(p,q) corresponds to the direc- 
tion of the point load and the second index refers to the component 
of the respective displacements and tractions. 

For the two dimensional isotropic plane strain case the 
fundamental tensors are [12] 


and 


' ■ ii(^l-vi r [If ^^ij '*,i'*,j' 


' "j^i * "i^j 


(3-17) 


The case of plane stress can be hardled through the use of an 
effective Poisson's ratio oiven by [v/(l+v)]. 

Due to the singular nature of Kelvin's problem at the point 
p(x) a small spherical region (circular region for the two- 


c - z 
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dimensional case) of radius e denoted surrounding the point p 
with surface S has been deleted (see Figure 3-1). A limit will 

C 

be taken by diminishing the radius e indefinitely. The external 
boundary of the region will be taken to be the surface S of the body. 



Figure 3-1. Develooment of Somigliana Identity 


Now Eouation (3-6) as apolied to the region R-R can be rewritten as 


f t.u^ dS 
S + S 


u^-t^ dS 


(3-18) 


S + S 


Noting that on the surface ^ ^ deleted 

region s * r; it is easily verified that the limit 


'Mm 

C-M) 


t.u^ dS « 0 


(3-19) 


The second integral in Equation (3-18) over may be written as 
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I u^t* dS » I [Il^(Q) - u.(p)jt*(p.Q)dS(0) 

Sg (3-20) 

♦ u^(p) I t*(p,Q)dS(Q) 

S 

e 

where point Q Is on the surface and point p(x) Is held fixed In 
the Integration. Since the displacements are continuous at p(x) 
the second Integral In Equation (3-20) has zero limit as e-K). And 

lim f u.t*dS * lim {u.(p) f t*dS} 
e-H3 ^ ’ c-o ’ ^ 

(3-21) 

- u,(p) SjjPj . Uj(p)ej 

Substitution of the above results Into Eouation (3-18) and taking 
the limits as e-*o 

I t.u* dS * Uj(p)®j + I (3-22) 

S S 

it ♦ 

Making the substitution that u. * 
taking each’ of the e. terms as independent the following vector 

V 

eouation is obtained 

Uj.(p) » f t^(Q)Uj^(p,0)dS(Q) - I u^.(Q)T^^(p.Q)dS(Q) (3-23) 

S S 

where u.(p) is the displacement vector at an arbitrary internal 
point p, and the vectors u^{Q) and t^(Q) are the displacement and 
traction vectors at a surface point Q. Equation (3-23) is the well 
known Somigliana identity for the interior displacement field 
generated by a full set of known boundary displacements and 
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boundary tractions (cf. Reference 31, p. 245). 

The interior stress state may be generated by differentiation 
of Equation (3*23) with respect to the load point p(x) and using the 
constitutive equation. Equation (3*2), and is given by 

<Tij.(p) = I t^(Q)D^^j.(p,0)dS(Q) - j u^(Q)S^.j(p,Q)dS(Q) (3-24) 

S S 

The third order tensor kernels and are given by [32] 
Djjij(p,Q) » (-“J Q^*2v)(5|^.r^j. + 5^.r^^ 




(3-25) 


^ ‘i/>k 

+ a v(n.r -r . + n.r .r . ) (3-26) 

• *J J 

- u-aviVij I 

In the above two equations o»l,2; 8=2,3; ■r*4,5 for two and three 


dimension problems respectively. The coima-differentiation is to be 

taken with respect to the coordinates of the point of integration, 

and since the stresses are found by taking derivatives of Equation 

(3-23) with respect to the solution point p(x), a change of sign 

arises from the identity (^ ■-“ ). The normal is evaluated at 

3x^ 3y,- 

Q(x). 

While Equation (3-23) reprecents the interior solution , its use 
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requires the simultaneous specification of all components of 
both the displacement vector and traction vector for all the 
boundary points. These six boundary conditions are twice as many as 
required by the uniqueness theoran. 

The essence of the boundary* Integral equation technique Is to 
allow the Internal point p(x) to pass to an arbitrary surface point 
PCj). yielding a set of Integral equations which can be solved for 
the unknown boundary tractions and displacements In a wel Imposed 
boundary value problem. It should be noted that the kernels of the 
Integrals in Equation (3-23) exhibit singularities as p(x) 
approaches the surface S, so that the resulting set of Integral 
equations will be singular. 

The limiting procedure for p(x) approaching a point P(x) on a 
smooth contour (two-dimensional case) Is Illustrated schematically in 
Figure 3-2a. This smooth contour permits the construction of a short, 
straight line centered at P(x) with length 2e. The Internal point 
p(x) Is placed at the boundary and the boundary Is augmented as shown 
In Figure 3- 2b. Point p(x) Is assumed to be at the center of the 
circle and afterwards the radius e Is reduced to zero. The point 


will then become a boundary point. The results of this limit for 
the Individual terms of Equation (3-23) are as follows 


11m u.(p) » u.(P) 

p(x) -► P(x) ^ ^ 


(3-27) 


11«n u.(Q)T..(p,Q)dS(Q) » 

p(x) - P(x) J ^ 


S 
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lim 

e-^0 


I “i' 


,T,,(1S + Hm 

s-s ^ 

£ 




For simplicity the second limit on the right hand side will be 
called I; ^ 

I . m j u,Tj,<lS . 11 J - j u,[|t |r^l 


♦ U-2vH|? n, - n,)] 


3r 


1 


dS 


j * 3x. "i'*'^ 4ff(l-v)r 
J 

knowing that ear and using a polar system of coordinates for 

simplicity,, the second term in the above limit will disappear as, 

3r ^ 3r ^ ar 3r ^ 3r 3r . « 

3x^ * 3Xj *'i 3X^ 3Xj 3Xj 3x^ ~ 

Since ■jj;' * 1 and dS * rde where the range on 9 Is -a^e^ir-a then. 


11m 

c-0 


I u,[(l-2v)«,_j ♦ 2 I 


This can be expanded and expressed in terms of 9 only, for instance 
when j = 1 


j ir-a 

I • < - I [uj(l-2v)+2UjCos^9 + 2u2Stn9cos8] 




9=-a 


.11m <- tu,(9)> 

- 9-0 I 2 1 j 


■ 


The same thing can be shown to apply for j » 2 The limit is then 
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11m u.(Q)T..(p.Q)dS(Q) 

p(x) - P(x) J ' 

' S 

, r (3-28) 

«-f Uj.(P) + ju.(Q)Tj..(P,Q)dS(Q) 

S 


where the Cauchy Principal Value of the above Integral Is taken. 
The llroix 


11m 

p(x) -*• P(x) 


jt^(Q)Uj.,.(p,Q)dS(Q) 

S 

[t-U.. 
j 1 JT 


lim 

e-0 


dS + 11m 
e-^ 




s-s. 


The second limit on the right hand side can easily be shown to be 
zero, hence 


lim 

pU) - P(x) 


t.(Q)U^.(p,Q)dS(Q) 


|t^-(Q)Uj^(P,Q)dS(Q) 


(3-29) 


Substituting these limtis. Equations (3-27), (3-28) and (3-29), 
yields the boundary constraint equation for the case of smooth 
boundary 


I 


u.(P) 



u^(Q)T^^(?,Q)dS(Q) 


» I t.(Q)Uj.(P.Q)dS(Q) 
S 


(3-30) 


The BIE for the three-dimen: onal case essentially has the same form 
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as Equation (3*30) and Its derivation may be found in reference [33]. 

For the case of the non smooth boundaries, the BIE, Equation 
(3-30), takes the form [34]. 

Cj,-(P)u^(P) + f u.(Q)Tj^.(P,Q)dS(Q) 


S 

» I ti(Q)Uj,.(P,Q)dS(Q) 


(3-31) 


S 


In Equation (3-31) c. . is a matrix of m^nbers and depends only up 
on local geanetry of the boundary at P(x). In practice, c^.^ is most 
expediently calculated by noting that an arbitrary rigid body trans- 
lation of R is a sol of Equation (3-31) [15]. If the body 
force, b., is considered. Equation (3-31) takes the form [23]. 

Cj^(P)u.(P) + I u.(Q)T^.^(P,Q)dS(Q) 


^ (3-32) 

« I t.(Q)Uj.(P,Q)dS(Q) + j b.(q)U^.(P,q)dV(q) 


S 


R 


Equation (3-31) can be viewed as the constraint equation re- 
lating surface tractions to surface displacements. In physical 
problems the tractions and displacements are not known concurrently 
over the entire surface. Thus the mechanism of solution is to regard 
Equation (3-31) as a set of coupled integral equations of varying 
types according as data appropriate to the traction, displacement or 
mixed boundary value problems are prescribed. 

Exact solutions to Equation (3-31) are, in general, not obtain- 
able and Equation (3-31) must be solved numerically for the unknown 
boundary data. Once this is done, the solution for displacement or 
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stress at any interior point p(x) of R may be obtained, in terms of 
the surface pair u^, t^, by comparatively simple quadrature using 
Equations (3-23} and (3-24) respectively. 

A closing discussion of this section is the treatment of 
Equation (3-32) where body force is included. The presence of body 
force b^(q) in this equation removes the problem from 'boundary only' 
category, assuming that the volume integral cannot be analytically 
performed. However, it has been shown in Reference [35] that it is 
possible to convert this volume integral to a surface integral when the 
body force, real or thermal equivalent, is derivable from a scalar 
potential whose Laplacian is at most a constant. Now if there exists 
a scalar t|i such that 

b. s “ . with ♦ ,.,* = k. (k^ * constant) 

1 ,1 ,n 0 0 

The volume integral in Equation (3-32) is convertible to an S 

integral and the BIE for this case may be written as [23] 


c..(P)u^(P) + [u^(Q)T-^.(P.Q) - lc.(Q) - ^(Q)n^(Q)} 

J 1 J J 

s’ Uj^(P,Q)]dS(Q) 

3r i 


(3-33) 


[t|;(Q) -^(P.Q) - (^*(Q)} r^j(P.Q) 


+ kQn^.(Q)r(P,Q)]dS(Q) 


where = ( 1 - 2 v)/ 16 t 7 u( 1 “v) 

Equation (3-33) is now amenable to a systematic 'boundary only* 
numerical treatment. 
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3.3 Formulation of the Numerical Solution Technique for 
Plane Problems 

As stated In the previous section* analytical solution to the 
BIE» Equation (3-31), are not generally available. A general scher • 
for obtaining numerical solutions for the case of plane elasticity 
is discussed in this section. This scheme is to replace the 
continuous definition of the boundary data in Equation (3-31) by 
some discrete variation that will reduce it to a set of linear 
algebraic equations. The boundary is approximated by a set of bound- 
ary segments over which u., t. vary in some specific manner. 

J M 

The approach followed here is the use of isoparametric shape 
functions for the representation of both geometry and boundary data 
[36]. The basis of the isoparametric method is the mapping of a 
planar curve to a standard interval by means of a fixed set of shape 
functions. The boundary is discretized into a total number M of 
boundary segments; each is defined by n nodes. The cartesian co- 
ordinates of each node are specified, and the cartesian 
co-ordinates of a non-nodal point of a segment are assuiued to be 
given by* (cf. Reference 36). 

X,(C) . N“(c) X? ,3.3,) 

in which N®(c) are the shape functions of intrinsic coordinate 
C(-l S 5 < 1) of the boundary element. 


* Generally for the isoparametric representation of geometry not 
higher than quadratic shape function are used [36]. 
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Boundary data can be modeled using a similar isoparametric 
representation as used for the geome:ry, Equation (3-34). The dis- 
placement vector of a point Q(x) on segment (o) is given by 


Uj(Q> 

UgCQ) 


N u 

- -s 


(3-35) 


where Ns 


rN'(5)n2(ON2(5) 

000 .... 

1 

o 

o 

o 

• 

• 

• 

• 

N^(c)N^(5)N^(C) .... 


and u^ is a vector contains the displacements of the nodal points of 
the segment; namely 


u"*" = {uJu?u? 
-a 111 


1 2 3 

U^U.Ua 

2 2 2 


} 


seg.a 


with similar expression for!tj(Q)| 


|t2(Q)l 


Writing Equation (3-31 ) for a nodal point P (see Figure 3-3) 

n 




M 

r 

a=l 


f u.(0)T..(P ,Q)dS(Q) 
; ' J' n 

S 


M 

Z 

0*1 


j t.(Q)Uj^(P^.O)dS(Q) 


(3-36) 


S 

0 

where the integrals of Equation (3-31) are now divided into the sum 
of M integrals along the surfaces S of the boundary segments. 
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Figure 3-3. Integration Along Boundary Segment {a) 


Based on the discretizing assumptions, the parametric 
functional representations for u^(Q) and t^(Q), Equation (3-34) are 
substituted into Equation (3-36) to obtain for a particular nodal 
point P 

n 



n I n 


M n ( \ 

+ r I 
0=1 0=1 ’ 


T..(P ,Q(0)N“(c)J(c)d5 

j * n 

s 

® (3-37) 


M 

r 


a»l 


(a, a) 


Uj^(P^.Q(0)N' 


(5) J (c)dC 


c 

"a 

where and are the displacements and tractions at 

local node o on the t.undary segment a. Suimation is implied on the 
repeated indices i and j. The range of the indices i and j is now 
{1,2}. j(c) is the Jacobian which more explicity is 
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Writing Equation (3-37) for every node and performing the inte- 
grations over each segment on S, including that segment which 
contains itself, a systan of 2N algebraic equations will result 
where N is the total nunter of nodes. 

Introducing the short notation 




(3-39) 


Sa 

Equation (3-37) can be written in the following form 
I " u ■ E e" t 

’ j=l o»l ^ na o»l c=l ’ no 

(n * 1.2 N) (3-40) 

The integrals in Equation (3-39) have the meaning that integration 
of the kernel function over the segment o for the nodal poi-nt n. 

Introducing a global numbering for the nodes along the boundary 
S, this leads to viewing the system of equations (3-40) in the 
following global matrix form 

cU + hU«6T (3-41) 

•w Ml ^ m W « 


or 

H y « G T (3-42) 

In Equation (3-42), the column matrices U and \ are of dimension 2N 
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a g 

and contain the elements u?» t?, 3 * 1,2 ,N, which are, 

respectively, the displacement and traction components at global 
node number 3. The square matrices H s c h and G are 2Nx2N 
elements In size and contain the integral coefficients of U and T 
as Indicated In Equation (3-40). These coefficients are sums of 

Integrals of kernel -shape function products, and generally evaluated 
using Gaussian quadrature formulae, c Is a diagonal matrix, which 
does not need to be determined exp 11 city. 

For the case when n*S , there Is no need to evaluate separately 
the Cauchy Principal Value of the Integral Equation (3-40), 

and the coefficients c- of the free term; so the coefficients In 
the H matrix can be calculated using the fact that the stress field 
corresponding to a rigid body translation is zero. Assuming unit 
rigid body translations (In the and X 2 directions). Equation 
(3-42) yields 

H I « 0 (3-43) 

where I Is a unit vector. Thus the sum of all the elements of H In 
any row ought to be zero, and the value of the coefficient on the 
diagonal Is simply equal to the sum of all the off-diagonal terms 
with sign changed. 

The numerical procedure for solving the BIE for plane elasti- 
city problems discussed thus far will now be applied to the case of 
assuming linear variation of the boundary functions. The linear 
boundary element scheme is simple and direct and offers a greater 
flexibility in solving general problems in solid mechanics. It is 
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quite adequate for practical applications when the 3IE method 
is used in conjunction with the experimental techniques for measuring 
in^lane displacement components. For these applications the region 
of interest is enclosed by a closed contour on which the displace- 
ment components at the nodal points are measured experimentally. 
Normally the region of interest is relatively a small area and, 
therefore, its contour may be approximated by a sequence of short 
straight line segronts. A c(xnputer program based on the linear 
functional variations along the boundary was used for the BIE calcu- 
lations for the example problems conducted in this dissertation. 

Good degree of accuracy has been achieved from using this procedure 
as will be shown in Chapter V. 

For linear elements (Figures 3-4), the variation of the 
functions u^(Q) and t^(Q) between nodes is assumed to be linear. 

The nodes are considered to be at the intersection between two 
straight line elements and marked as extrene nodes. There are two 
shape functions N^ = N^(£) • • y (^-1) 

and = N^(S) « j (S+D 
1 

and the Jacobian is simply J(f) * where 1^ is the element length. 



Figure 3-4. Linear Boundary Element 
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The assembled equation for node n» Equation (3-37), nuy be 
written in a matrix form for the case at hand as 


s''!;'’ * 

(a • 1,2 M) 

In Equation (3*44) a is a (2x4) submatrix defined by 


(3-44) 


5-10 


sj 

."Vi . 

L S< 

' N^T^dC i 
1 Sc 

nH-.,ds i 

* ' J 

i Sc 

■ N^TjjdS 

J 

-si 

"Vi j 

1 ! S< 

f"V'i , 

T 1 S< 

1 ' 

N%-dc 1 

C£. , J 

T 1 Sc 

N^T22d« 

J 


with a similar definition for B 

-no 


The elenents of the submatrix 


a^^ are integrals along the element o corresponding to node n. 

Note that N^ and are functions of s and T. . s T..(P ,Q(c)); 

I J I J I 

* the Jacobian J( 5 ) and is the length of segment. 

It is more convenient to rewrite Equation (3-44) n terms of the 
nodal displaceiT«nt and traction vectors in the following forms 

-np - -HP - (3.45) 

(p * 1.2 N) 

where N is the total number of nodes « the total nunter of segments, 
and the vector u^ is the displacement vector at node (p); namely 



with a similar definition for t*^. The coefficient A in Equation 

• -no 

(r-45) is now a (2x2) submatrix given by 
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which contains integral over the two adjacent segments meeting at 
node p. When p«I, N is used instead of (p*l) in the above express- 
ion. 


At any node there will be a single unique vector for the nodal 
displacement; the same thing is not always true for the traction 
vector at the common nodal point. To properly account for discon- 
tinuities of nodal tractions a zero-length segment is to be 
introduced at that nodal point. Riccardella [37] has shown for the 
case of linear functional variations that zero- length segments can be 
added at any point where a step change in boundary conditions occurs, 
and ze* > values can be assigned to the integrals associated with 
these zero- length segments. 

Writing Equation (3-45) for every node, a system of algebraic 
equations will result, 

(3-46) 

or 

[A, 1 {u“> • [3„ J (t®} (3-47) 

• rip • -HP 

where 

A , “ A when n f o 

A “A + c’^ when n “ 

-np -HP 


0 


(3-48) 
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The Integrals for [A ] and [B ] may be evaluated in closed form. 

•“HP "•HP 

Again there is no need to calculate separately the diagonal coef> 
ficients of the matrix [A ] as discussed earlier in this section. 

-HP 

Assuming a unit rigid body translation in both and x^ direction, 
i.e., {jj®} = I, the corresponding boundary traction is {t*’} s 0, 
then from Equation (3-47) the value of the diagonal coefficients is 



N 

- l 
p-1 



(3-49) 


Pfn 

Any well-posed boundary value problem will have 2N knowns and 
2N unknowns. Once the knowns of the boundary data have been 
specified, the system of equations. Equation (3-42), can be re- 
ordered in such a way that all the unknowns are on the left side, i.e.. 


la 


= Y 


(3-50) 


The vector X will contain 2N unknown nodal values of traction or 
displacement. The vector Y contains the product of the calculated 
integrals and the corresponding known boundary data. The system of 
equations (3-50) is to be formed directly, and thus the storage of 
matrices H and G has been avoided. Only the matrix of coefficients 
of the unknowns at the nodes need be formed and stored in a square 
array (matrix K). The elements of the other matrix of coefficients 
to be multiplied by known quantities may be inmediately multiplied 
by such quantities as they are generated and stored in a column 
vector and not as a separate square array. 

The order of magnitude of the terms in matrix H is approximately 
(u) times that of the terms in matrix 6 (where y is the shear modulus). 
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Thus during the formation of matrix K, the coefficients of the un- 
known nodal tractions have been scaled using the shear modulus to 
maintain the same order of magnitude for tne unknowns. The systen 
of Equations (3-50) is numerically st^le and is now directly solv- 
able by standard reduction schane (Gaussian elimination). 

While one of the major potential numerical advantages of in- 
tegral equation methods is reduced problem dimension by solving 
boundary equations, another major feature is the internal solution 
capability. Once the system of Equations (3-50) is solved, the now 
completely known discretized boundary data may be used to evaluate 
the displacements and stress tensor at any interior point p(x) by 
direct integration of the identities (3-23) and (3-24) with the 
boundary- integrals discretized in the same manner as Equation (3-37); 


u.(p) * z r f U...(p,Q(i))N®(c)0(5)dC 

^ C»1 0*1 ’ J 


- z " f T.,(p.Q(5))N®(c)J(UdC 

a»l 0*1 ’ J 


..(p) • Z Z t.^®*®' [ D. .,(p,Q(5))N®(c)J(5)dC 


a-1 ' 

, C*1 0*1 


where suimiation is implied on repeated i or k indices. Again the 
indicated integrations may be evaluated numerically (using Gaussian 
quadrature formula) and u.(p) and c..(p) are found for any interior 

W * V 
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point. By keeping the discretization on the boundary very high 
numerical resolution may be obtained, from Equations (3-51) and 
(3-52), a feature quite unlike finite difference and finite elenent 
analyses. When high resolution is required or only a limited 
number of internal point solutions are desired the integral equation 
method is significanlty faster than finite element solutions. 

It should be mentioned here that the calculated values of the 
stress tensor at interior points very near to the boundary are not 
very accurate, since the boundary is defined only by the coordinates 
of the nodes, and traction and displacement data are av:^ilable only 
at the nodes. It has been found that an interior point can Le with- 
in a distance equal to the length of the nearest boundary segment 
and still maintain ordinary accuracy, i.e., accuracy comparable to 
that available for the other interior points at distance from the 
boundary more than the length of nearest segment. 

Not being able to take p(x) arbitrarily close to the boundary is 
not a limitation, however, since stresses right on the boundary may 
be obtained with no reference to the interior field. The boundary 


stress field may be obtained directly from knowledge of both 


tractions and boundary displacements. Since at this stage u. on S 

du. 

is known at the nodes such that can be obtained mmierically 

(and sometimes experimentally). Noting further that 


du^ 
ds p 




* “i,2"l 



(3-53) 


where q is the unit tangent to S at P, and the tangential derivative 
du^ 

is assumed keeping the material on the left while moving in a 
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direction of positive S. Equation (3*53) together with Hooke's law. 
Equation (3-2), and the boundary traction t^ » ®ij”j ^ consti- 
tute sufficient information to determine the state of stress, o. ., 

* J 

at each node. The above relation can be casted in the following 
matrix form [18] 

BJ = M (3-54) 

in which 



The square matrix B contains the elastic constants and the components 
of n^ and q^. After obtaining numerically (or sometimes experimen- 
tally) the components u, and u- at a given node point P, the 
systan of Equations (3-54) becomes ready for convenient computation. 
The derivatives u. plus t. at P complete the column M such that 
the stress components o.. on S are obtained as the solution of 

* w 

Equation (3-54) for each node. Along with the stresses, the dis- 
placement gradients u. . are also obtained, although the later are 
not normally of interest. 

The BIE solution method can easily be used to study problems 
with more than one surface, such as the case of a body with holes 
illustrated in Figure 3-5. In order to define an external or 




53 


internal boundary, direction of the normal needs to be identified. 
This can be done for two dimensional problems by defining a number- 
ing scheme as the boundary is traversed in the counterclockwise 
direction keeping the body material to the left. 

Finally, a basic conceptual advantage of the BIE method is 
that it applies to three dimensions as well as two with similarly 
defined integral equations arising in each case [13]. The proce- 
dure for three dimensional problems is similar, except the number 
of equations and unknowns is now larger, and evaluation of the 
coefficients involves surface integrals over the surface segments 
rather than line integrals over line segments. 



54 


IV. THE EXPERIMENTAL PROGRAM 
AND DATA ACQUISITION 

4.1 Experimental Equipment 

The experimental arrangement for making a double exposure 
specklegram is shown in Figure 4«1. It is, in fact, a simple 
camera configuration. Monochromatic and coherent laser light is 
used to illuminate the specimen which provides for a speckle 
appearance on a diffuse surface. Light reflected from vne object is 
imaged onto high resolution photographic film. Thin emulsion holo- 
graphic films such as Agfa-Gevaert holotest 10E56 are commonly used. 
For the point-by-point data analysis, a set up similar to the one 
shown in Figure 2-6 may be used. A helium-neon laser with one to 
fifteen milli-watts of output will provide adequate illumination. 

To obtain good resolution of changes in displacement with respect to 
location, a small diameter input beam is preferred. The diffraction 
halo with fringes is observed on a screen. The procedure for deter- 
ming the in.pl ane displacement components at a point was discussed 
in Section 2.3 of Chapter II. 

A computer aided data reduction system was developed primarily 
for apoli cation in speckle interferometry data analysis [5]. This 
system can be used as a mean for the displacement data acquisition 
for the example problems discussed in this dissertation. The basic 
idea of the operation of this system is discussed in the next section. 
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4.2 Computer«Aided Data Reduction System 

The basic system consists of two par* which are an optical 
recording system and the computer hardware used in the numerical 
analysis as shown in Figure 4-2. 

A typical fringe pattern used for analysis is shown in Figure 
4-3. This type of fringe pattern is produced when a laser beam is 
passed through a photographic transparence as discussed earlier. 
Interference fringes are a measure of surface deformation of a body. 
Be'*ause only a small area is illiminated by the laser, the infor- 
mation in the data analysis yields the displacement at a print. 
Therefore, to obtain a complete map of the surface, each point on 
the film will have to be illuminated. This is accomplished by an 
x-y translation stage in which the transparency has been mounted. 

The stage has the capacity to translate 6 inches in each direction in 
0.001 inch increments with a .0001 inch repositioning accuracy. 
Synchronous stepping motors provide he control for the translation 
directions. A photographic plate is mounted in the window of the 
translation stage and* a He-Ne laser beam passes through a point on 
the film and the film is translated relative to the sf^tionary laser 
beam. 

Location of a point in the window is controlled by the stepping 
motors which are in turn controlled by the computer. Thus, the 
computer can control the *ranslation distance with a iTiinimum inter- 
val of 0.001 inch. 

••t that is passed through the film then forms a diffraction 
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Diffuser 



Storaoe 

Display 


Figure 4-2. Computer Aided Data Analysis System 
for L..S Speck U Interferometry 






Figure 4*t 3. Typical Speckle Photography Fringe 


Pattern Used In the Optical Data 
Analysis System 
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pattern which Is imaged by a vidicon camera. The Image on the 
camera can then be viewed on a TV nranitor or stored in memory in 
the ccanputer. 

A POP 11/40 miniconputer is used to control the motion of the 
translation table and process the video signal. Operation of the 
system begins with a coninand from the computer to move the table to 
a specified coordinate position in the table window. With the table 
in a specified position, the amiputer sends a conanand through the 
video interface to store the vidicon imager in the scan conversion 
meiTOry. A ccnmand from the conouter t^en converts the information 
stored in the scan conversion memory into digital form. The 
complete image is not processed but only selected lines of video 
Information are processed because of the requirements in the data 
analysis. This process completes the data acquisition at a point on 
the film and the computer then locates a new coordinate position and 
the process is repeated. 

Data analysis consists of determining the locus of points of 
minimum intensity which is referred to as fringe. However, this 
analysis is complicated by the fact that the film transparency is 
illuminated by a coherent lighc source which produces laser speckle. 
The presence of this speckle is to superin^ose a noise-like signal 
on the spatial information which degrades the image. A speckle 
averaging technique is used which reduces the effect of noise in 
locating points of minimum intensity. Fringe separation is deter- 
mined by first locating the point of maximum intensity from the 


original PAQ6 *8 
OF quality 
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scan average. Then the first minimum point is located to r;.ch side 
of this maximum and the difference in location determines the fringe 
spacing. 

Several examples will now be considered vtd valt^s of the 
measured displacenents will be presented. The measured values of 
displacanents are used as data input to numerically calculate 
stresses and displac^nts through ^e BIE method as will be dis> 
cussed later in Chapter V. 

4.3 Example 1 - Flexure of a Cantilever Beam With a Transverse 
£nd Load 

A beam with tite end x « 5.0 inches fixed and the end x « 0 has 
an end load directed along the negative y-axis was considered and 
shown in Figure 4-4. Material for the beam was 1/4 inch plexiglas 
witti £ * 4.56x10^ psi and v « 0.37. Fixed end conditions for the 
beam were obtained by bonding one end of the plexiglas beam to an 
alixninum block. The transverse end load was governed by the 
deflection of the free end (because of the constraints on the pcint- 
wise filtering of the experimental data). The end load was applied 
with a micrometer screw and the deflection was recorded with a dial 
indicator. Free end deflection was 5 * 0.005 inches which corre- 
sponds to an end load Vl^ * 1.14 lb. 

Displacements along the free end were measured experimentally, 
using the laser speckle technique, and compared with the theoretical 
values predicted by the stress function solution. The comparison 
as depicted Figures 4-5 and 4-6 shows a favorable agreement. 

Also experimental and theoretical values of the displacement along 
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the neutral axis of the beam were 1n a very good agreement as shown 
previously in Figure 2-8. These favorable comparisons between 
experiments and theory, at the free end and along the neutral axis 
of the beam, encourage the use of such experimental technique for 
measuring the in.pl ane displacement components along any specified 
contour. 

For the region of interest shown in Figure 4-4, displacanent 
components were measured along the fixed end (x * 5.0 inches) which 
are zero and along the section x * 3.0 inches. Data analysis for 
this line yielded the experimental values of displacement 
components shown in Table 4>1. This completes the displacement 
data needed for the BIE solution for this region. 

Also the displacement components at some chosen points inside 
the region of interest along the x-axis were measured and shown 
in Table 4-2. These values will be used later for the comparison 
with the BIE calculations as will be discussed in the next chapter. 
It should be mentioned here that for points at distances greater 
than 3.6 inches from the free end fringes could ndt be observed 
because of the smallness of the displacements at these points. 
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TABLE 4-1. EXPERIMENTAL VALUES OF THE DISPLACEMENT 
COMPONENTS ALONG SECTION x * 3.0 INCHES 
FOR THE CANTILEVER BEAM EXAMPLE 


y 

(inches) 

uxlO^ 

(inches) 

vxlO^ 

(inches) 

-0.5 

-4.633 

-10.916 

-0.4 

-3.531 

-10.868 

-0.3 

-2.715 

-10.890 

-0.2 

-1.666 

-10.522 

• 

o 

1 

-0.913 

-10.435 

0.0 

0.000 

-10.475 

0.1 

1.095 

-10.418 

0.2 

1.788 

-10.688 

0.3 

2.437 

-10.559 

0.4 

3.500 

-11.100 

0.5 

4.633 

-10.916 


TABLE 4-2. MEASURED DISPLACEMENTS ALONG THE X-AXIS 
FOR THE CANTILEVER BEAM EXAMPLE 


X 

u 

o 

X 

> 

(inches) 


(inches) 

3.2 

0.0 

-8.469 

3.4 

0.0 

-6.686 

3.6 

0.0 

-5.326 
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4.4 Example 2 - Circular Disk Under Diametral Compression 

For this second example problem a circular disk of 3 inches 

diameter made of 1/4 inch plexiglas was loaded in compression by 

two diametrically opposing nearly concentrated forces located on 

the top and bottom of the disk as shown in Figure 4-7. Displacement 

components at the nodal points along the contour of the region of 

interest, shown In Figure 4-7, were measured experimentally. The 

vertical components of Idie displacement along the horizontal dia-. 

meter were constant. A zero value was assigned to these ccr.ponents 

because of the symmetry of the problem about the horizontal diameter. 

The horizontal components of the displacement along the vertical 

diameter were measured and are zero. Data analysis for line y=d/4 

yielded the experimental values of displacement components along 

that line which are plotted in dimensionless form in Figure 4-8. 

The horizontal displacement component at node 11 and the vertical 

displacement component at node 26 were measured to be u^^ « 

-5 -5 

19.7x10 and ® -36.9x10 inches respectively. Data input for 
the BIE solution on prescribed nodes along line y=d/4 were taken 
from the best curve fit of the experimental data shewn in Figure 4-8. 

Both the horizontal displacement components along the horizontal 
diameter and the vertical displacement components along the 
vertical diameter were also measured and will be compared with the 
BIE calculations as will be shown later in Chapter V. 
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4.5 Examole 3 - A Plate with a Central Circular Hole 
sub.iected to linifonn Axl^l lensly 

The geometry of the 0.228 inch thickness plexiglas plate used 
to study this classic stress concentration problem Is shown in Figure 
* . The ratio of hole diameter to plate width is 1/9. A load P of 

.lagnitude 140 pounds was applied in the y direction through 13 
equally spaced small holes to assure uniform tension at sections of 
the plate far from the loading ends. 

The region of interest selected for the stress analysis using 
the BIE method is the portion of the plate within a concentric 
circle of a alameter four times the hole diameter. Due to symmetry 
of the problem, only one quadrant of this region is considered. 
Displacement components at 19 points, of equal angular intervals, 
along the circular contour of the region of interest (nodes 11-29) 
were measured experimentally. Table 4-3 is a display of these 
measured values i.i dimensionless form. Again, a zero value was 
assigned to the measured constant value of the vertical displacement 
components along the x-axis, which is used as a reference for 


measurement of the vertical displacement components. 




Figure 4-9. A Plate with a Hole 
under Axial Loading 
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TABLE 4-3. MEASURED DISPLACEMENTS ALONG 
THE CIRCULAR CONTOUR OF THE 
REGION OF INTEREST FOR EXAMPLE 3 


Node 

“/“u 

v/*29 

11 

1.000 

0.000 

12 

0.992 

0.075 

13 

0.975 

0.150 

14 

C.941 

0.224 

15 

0.899 

0.299 

16 

0.846 

0.373 

17 

0.737 

0.446 

18 

0.719 

0.518 

19 

0.652 

0.589 

20 

0.579 

0.656 

21 

0.506 

0.720 

22 1 

0.433 

0.781 

23 

0.362 

0.834 

24 

0.298 

0.882 

25 

0.233 

0.924 

?t 

0.171 

0.956 

27 

0.112 

0.980 

28 

0.056 

0.994 

29 

0.000 

1.000 


-5 

Measured values of » -14.2x10 inches 
and 

V 2 g • 35.7x10“^ inches 
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V, RESULTS AND DISCUSSION 

5.1 The Coupling of the BIE anj Laser Speckle Techniques 

Experimental displacement measurements coupled with the BIE method 

represent a useful solution procedure which extends the capabilities 
of coherent cities measuranents in calculation of stresses or 
strains in engineering problems. Data analysis in speckle inter- 
ferometry yields a map of surface displacements along the contour 
of a region of interest. These experimental values are used as 
data input to numerically calculate, with a hi^ degree of 
resolution, strains and stresses at any desired location inside 
the region of interest through the BIE method. 

In order to demonstrate the coupling of. the two solution 
techniques as applies to linear elasticity problems, several examples 
are presented in the following sections. The nimierical results arr 
obtained through the use of a BIE computer program based on linear 
functional variations along the boundary segments. Displacement 
data input to the computer program are the experimental values of 
the example problems discussed previously in Chapter IV. 

5.2 The Can ci lever Beam Problem 

The cantilever beam problem was analyzed using the BIE model 
for r‘'e region of interest shown in Figure 5-1. Sixty straight 
line segments of equal length were used. Nodal points 1-21 are 
specified to be traction free, points 21-31 were measured dis- 
placements which are zero. Points 31-51 are traction free and 
points 51-1 are measured displacements uti’-zing the laser speckle 




Figure ?>-l. DIE Model for the Region of Interest 
of the Cantilever Beam Problem 


W 
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technique of pointwise filtering. These measured values were 
listed previ>usly <n ladle 4-:. Especial attention was given in 
the nunerical procedure to the tractiCT discontinuities at the 
comer points as discussed previously in Chapter III. 

Table 5-1 is a display of the calculated BIE results for the 
stress components at sections x«3.0 and 5.0 inches. The results at 
section x*3 inches were compared to the stress function solution 
and to the bending monent for a load » 1.14 Ibf. Integration of 

the stress components over this section yielded a bending moment 
deviation of -0.58S. The same results were compared at the fixed 
end (x>5 inches). At this section the bending moment and shear 
force dev* * from the static values by 6.49« and 7 6 % respectively. 

As antici;^ the stress distribution at the fixed end section 
is different t *ra the stress function solution. Since exact 
boundary conditions were specified for the BIE solution, it is 
logical to theorize that the results obtained from the combined 
laser speckle displacement measurements and BIE solution for this 
problem are more reliable than the stress function solution. 

Internal displacements at three points along the neutral axis 
of the beam were also computed and compared to the n^asured values. 
The comparison is shown in Table 5-2. 

The results of the combined laser speckle displacement measure- 
ments and BIE solution procedure compared favorably with the 
theoretical and the experimental results and static equivalent 
check. These results are particularly good in view of the small 



TABLE 5-1. RESULTS OF THE BIE CALCULATIONS 
FOR THE CANTILEVER BEAM PROBLEM 


y 

(Inches) 

0^ (psi) 
at x«3.0 Inches 

(psi) 

at x»5.0 Inches 

Tjjy (psi) 

at 

X » 5.0 in. 

Bir 

Theoretical 

■LUfli 


MEM 

76.04 

73.87 

127.03 

123.12 

14.61 


57.75 

57.46 

83.31 

85.76 

16.15 


41.34 

41.04 

57.80 

68.40 

-0.39 

IH 

24.60 

24.62 

34.48 

41.04 

-3.27 


8.09 

8.21 

11.44 

13.68 

-5.04 

BEH 

-8.09 

-8.21 

-11.44 

-13.68 

-5.04 

-0.15 

-24.60 

-24.62 

-34.48 

-41.01 

-3.27 

-0.25 

-41.34 

-41.04 

-57.80 

-68.40 

-0.39 

-0.35 

-57.75 

-57.46 

-83.31 

-85.76 

15.15 

-0.45 

-76.04 

-73.87 

-127.03 

-123.12 

14.61 


tn 
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applied transverse load W^. 


Table 5-2. COMPARISON BETWEEN CALCULATED 
AND MEASURED VALUES OF THE 
BEAM DEFLECTION ALONG THE X-AXIS 


X 

(inches) 

V X 10^ 

(inches) 

BIE 

Measured 

3.2 

-8.57 

-8.47 

3.4 

-6.89 

-6.69 

3.6 

-5.37 

-5.33 


The ability to partition a section of a geometry has an 
additional advantage in experimental data displacsient measure- 
ments. For e;ample, in this sample problan, the free end load was 
governed by the displacement measurenents of the free end. The 
ability to partition the geometry will allow displacements within 
■epginnc tn pxrpod fhp mpatnrahTp data if Only information is 
obtained in regions of high stress gradients. 
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5.3 The Circular Disk Problem 

This problem was studied using the BIE model Illustrated In 
Figure 5-2. The boundary of the region of Interest Is represented 
by thirty straight line segments. Due to the symmetric nature of 
the problem, boundary conditions for nodal points 1-11 are 
specified as zero traction In the x-dl recti on and zero displace- 
ment In the y-dl recti on. Furthermore, nodal points 26-1 are 
specified as traction free In the y-dl recti on and zero displace- 
ment In the x-di recti on. Points 11-17 are traction free and points 
17-26 are specified displacements. Displacement Input data for 
the BIE solution were obtained from the best curve fit of the 
experimental data previously shown In Figure 4-8. 

The numerical results for the normalized stress distribution 
along the horizontal diameter of the disk were compared to the 
theoretical results of Reference [38]. The comparison is very 
good as depicted In Table 5-3. Also, the results of the calculated 
displacements along the horizontal and vertical diameters compare 
favorably with the experimental results as shown in Figures 5-3 
and 5-4 respectively. 



Figure 5-2. DIE Model for the Region of Intere«:t 
of the Circular Disk Problem 
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TABLE 5-3. COMPARISON OF THE BIE CALCULATIONS AND 

THEORETICAL SOLUTION ALONG THE HORIZONTAL 
DIAMETER FOR THE CIRCULAR DISK PROBLEM 


Dinensionless 

Distance 

X 

Diirensionless Stress ^ 

m 

d7I 

BIE 

Theoretical 

0.0 

-1.941 

-1.911 

0.1 

-1.890 

-1.860 

0.2 

-1.743 

-1.719 

0.3 

-1.524 

-1.506 

0.4 

-1.266 

-1.257 

0.5 

-0.996 

-0.993 

0.6 

-0.738 

-0.740 

0.7 

-0.501 



0.8 

-0.299 

1 

-0.309 

0.9 

-0.126 1 

-0.141 




0 2 0.3 0.4 


Figure 5- 


Ditnensionless Distance 


. Horizontal Displacement Components 
Along Line y * 0 




Uimensionless Ulsplaceinent (v/v 
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Figure 5-4. Vertical Displacement Components 
Along Line x « 0 
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5.4 The Stress Concentration Problem 

The example problem of the plate with a circular hole under 
uniform tension was selected to demonstrate the application of the 
coupled technique of laser speckle interferometry and BIE to solve 
a typical stress concentration problem. 

Figure 5*5 illustrates the BIE idealization for the symmetric 
portion of the region of interest used for the stress analysis of 
the problem. Twenty-seven straight line segments with twenty-nine 
nodes were used. Nodes 1-10, along the hole surface, are specified 
to be traction free. Oisplacanent components of nodes 11-29 are 
measured values as listed previously in Table 4-3. 

Internal stresses were computed at a series of points along 
the x-axis. The numerical results, normalized by the resultant 
applied load divided by the net cross-sectional area, were compared 
to the available closed form solution for a similar problem of an 
infinite plate as shown in Table 5-4. It is of most interest 
to notice that the stress values predicted by the numerical 
solution are larger than the corresponding values predicted by the 
closed form solution. However, this is expected since the 
numerical results are for a plate with finite width. The boundary 
stress at node 10 was also calculated and the computed elastic stress 
concentration factor, maximum calculated stress divided by the 
nominal stress, was 2.66 as compared to 2.67 given In 
Reference [39] for the plate with finite width. 



ide 

0 


Figure 5 


Twenty-Seven Segment BIE Model 
Used for the Stress Concentration 
Problem 
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TABLE 5-4. COMPARISON OF THE BIE CALCULATIONS AND 
THEORETICAL SOLUTION ALONG THE X-AXIS 
FOR THE STRESS CONCENTRATION PROBLEM 


Dimensionless 

Distance 

x-R 

R 

Stress Concentration 

BIE 

Theoretical 

i 0.02 

2.63 

2.55 

1 0.03 

2.55 

2.49 

0.04 

2.49 

2.44 

0.05 

2.43 

2.39 

0.08 

2.23 

2.25 

0.10 

2.20 

2.17 

0.20 

1.37 

1.84 

0.40 

1.49 

1.46 1 

0.60 

1.29 

1.27 ! 

0.80 

j 

1.17 

1.15 

1.00 

1.10 

1.08 

1.50 

1.01 

0.99 


Oi 


R is the radius of the hole 
Nominal Stress * 204.8 psi 



85 


VI. CONCLUSIONS 

The general theory for laser speckle shearing interferometry 
which includes both the contributions of in-plane and out of plane 
deformations was developed. Laser speckle interferometry provides 
-a ti. que to directly measure in.plane displacement components 
without the need for vibration isolation. This experimental technique 
was coupled with a numerical boundary- integral solution technique 
to ca'.culate stresses and strains at any point for a deformed 
solid utilizing the measured boundary displacement data. 

The analytical formulation of the boundary- integral equation 
method for two and three dimensional elasticity was reviewed in 
detail. The linear boundary element scheme used in the numerical 
formulation for plane problems is simple and offers a great 
flexibility and ease in solving many important solid mechanics 
problems. 

The numerical results of the sample problems presented 
illustrate the successful coupling of laser speckle interferometry 
and boundary-integral equation techniques as applied to elasticity 
problems. High accuracy was attainable by using a relatively 
coarse level of boundary discretization. The ability to partition 
a section of a geometry has an additional advantage in experimental 
displacement measurements. It may allow displacements within 
regions rather than the particular region of interest to exceed the 
measurable data. The coupling of laser speckle interferometry and 
the numerical boundary- integral technique becomes more compatible 
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when the displacement data can be stored in the computer directly 
without the need of time consairi-n^ dna'anelysls. The coupTed 
technique may be applied to axially syrnnetric and fracture mechanics 
problems. In particular, fracture mechanics applications appear 
to be a well posed problem using this coupled technique. A closed 
contour surrounding the crack will form the measured boundary 
displacements and the boundary-integral solution technique used to 
calculate internal values. Thus the fracture mechanics problems 
form the basis for future research work. 
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